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ABSTRACT 


The effect of an ice-cover on transverse spreading of a tra- 
cer released at a constant and steady rate into a turbulent channel 
flow is investigated. Velocity and tracer concentration measurements 
are conducted at various sections downstream of the tracer source in 
both a straight and a meandering flume for open-water and ice-cover 


conditions. 


Transverse diffusion coefficients (E,) in the straight chan- 
nel tests are computed by the method of moments and by a formula based 
on similarity considerations. Transverse mixing coefficients (Ey in 
the meandering channel tests are evaluated using the moments and inte- 
gral methods based on the convective-diffusion equation expressed in a 


meandering coordinate system. 


It is shown that the transverse diffusion coefficient for 
ice-cover conditions in a straight channel is approximately 50% smaller 
than the corresponding coefficient for open water conditions but when 
these coefficients are normalised with the product of shear velocity, 


u, , and hydraulic radius, R , they do not differ significantly. 


* 


Analysis of the available data on transverse diffusion does not indicate 


a clear dependence of E,/u,R on W/R , where W is channel width. 


The magnitude of the transverse mixing coefficient in meander- 


ing channel tests and for ice-cover conditions is found to be approxi- 
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mately 4 times smaller than the corresponding coefficient for open- 
water condition. Values of E wg! Uae are higher for open-water than 
for ice-cover conditions for similar flow conditions. The longitudinal 
variation of the transverse mixing coefficient is found to be strongly 


related to the growth, decay and reversal of the spiral motions. 


Field investigations of transverse mixing in the Lesser Slave © 


River confirm the reduction in mixing capacity due to an ice-cover. 
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CHAPTER 1 


INTRODUCTION 


1.1 General 


In recent years the ecological balance of many rivers has been 
seriously upset by unregulated disposal of pollutants. Because rivers 
perform functions that are essential to man's survival pollution control 
presents itself as a major technological challenge. Dilution of a ple 
lutant by a river is rarely a satisfactory solution for a major a ee 
tion problem, yet most solutions involve dilution processes. Continued 
study and research on mixing of pollutants in natural rivers are, there- 


fore, important. 


Although a good amount of knowledge on mixing in open channel 
flows has been accumulated in the past two decades, there is still room 
for enlarging our understanding of the processes involved and particu- 


larly for more reliable predictions of the spreading of pollutants. 


Transverse (across channel) mixing is important in rivers 
because it exerts a large influence on the rate of longitudinal disper- 
sion (Fischer, 1967). Knowledge of the transverse exchange (diffusion 
or mixing) coefficient is also a prerequisite in estimating the crossing 
distance and the distance required to achieve a given degree of mixing. 
Considerable practical importance therefore attends the estimation of 


rates of transverse mixing of pollutants especially in areas where 
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location of water intakes is crucial. 


Most rivers in Canada are ice-covered for nearly 5 months of 
the year and pollution problems in these rivers are almost always most 
critical during winter due to generally low winter flows and suppres- 
sion of oxygen uptake from the atmosphere. However, no laboratory or 
field measurements of the transverse mixing characteristics of ice- 


covered channels appear to have been reported so far. 


a2. Objectives 


It is the principal objective of this investigation to pro- 
vide information on transverse diffusion and mixing in straight and 
meandering ice-covered channel flows. The specific objectives are: 


(i) to develop working equations for evaluating the transverse 
diffusion and mixing coefficients in channels of arbitrary 
shape, 


(ii) to obtain experimental data in the laboratory and field 
for corresponding open and ice-covered channel flows, and 


(iii) to analyse the experimental results and evaluate the effect 
of ice-covers on the diffusion and mixing phenomena. 


Experiments were, therefore, designed to fulfil these objectives. 


1.3 Structure of Report 


This section deals with the form of presentation. 


In Chapter 2 previous studies on transverse diffusion and 


mixing are reviewed and simplifications of the 3-dimensional convective 
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diffusion equation into forms amenable to solution and the limitations 


of some known solutions are discussed. 


Chapter 3 deals with the analytical considerations. Two 
methods are proposed for evaluating the transverse mixing coefficient 


in meandering and straight channels of arbitrary cross-section. 


Chapter 4 summarises the experimental arrangement and proce- 
dures while the experimental results relating to both straight and 


meandering flumes are presented and discussed in Chapters 5 and 6. 


Chapter 7 contains the experimental procedure, results and a 
discussion of these results for field tests on the Lesser Slave River 


(Alberta). 


Overall summary and conclusions of both laboratory and field 


investigations form the subject of Chapter 8. 


enobsetiot( sit hes satsvice of eldgnemeemioy oad: 


sbseeuoerh era rokefce sont 


owl .eackissehbienos ci ans daiv atest aa. 


taotaitisos gnivim seravansrs PP gaideuievs 108 | 


.roligae-acoxo yrersidys to atone hianerte bers 


7a oe: 


y ’ A 1s i Pp bie 
“soo Tq bas Jnesmagnn IG fstqani 19¢x9 eile) agul relnnila z sass : lee 
bas sigiexte ified of grttokon cates ShLTLTiY vit ‘sini 


= 
-0 bes ¢ etstqst) ried baaervak tb Hees betnenaed sis: Joanie an 


e t 
s bos #32uee1 .siubacorq ingnsmiyegqxs sd entsones,§ tadqudo 


avi sual 


y 4 fer] 7 - e ¢ \ 
2 sstast $3 no e@taax blak? 202 adiveor bears 36 ie 


Lisl bas yroserodsl Haod 26 eiokaulanes td vtemme, iad 
NS 
8 eetqeay) io Jsatdue ada ee 


- : Z| 


CHAPTER 2 


REVIEW OF PAST STUDIES 


2.1 The Convective-Diffusion Equation 


The study of turbulent mixing involves the analysis of mat- 
erial transport by a turbulent fluid. A comprehensive approach to the 
problem demands both careful experimentation and perp hse: theoreti- 
peak ne iclite The usual theoretical approach is by way of a mass bal- 
ance or conservation equation for the quantity being transported. The 
principle of mass balance for some scalar quantity, c , (herein called 


a tracer) in an arbitrary and stationary control volume states: 


Rate of mass increase within the control volume = 
rate of mass inflow into the control volume + 


rate of mass production with the control volume 


If it is assumed that the tracer undergoes only the process 
of dilution and that its fluid properties are identical to those of the 


ambient fluid then the principle of mass balance may be stated mathe- 


matically as follows: 


P) 3 a oe d (pc) 
ory (ogee Ox, (ou, c) = x, om ax, ) noe1) 


i 
where Xs is the co-ordinate in the i" direction: i.= 1, 2,33 


u is the instantaneous fluid velocity in the x, direction; po is 


i 
the fluid density; De is the molecular diffusion coefficient; t is 
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time; and c is the tracer concentration expressed herein as a mass 
ratio. Equation 2.1 describes the movement of the tracer completely. 
It is valid for turbulent flows and provided it can be solved would 


yield the instantaneous values of c. 


In general, however, the fluctuations of u, in turbulent 
flow are not known and in many cases knowledge of the fluctuations of 
c is not required. For turbulent flows, it has therefore become 
customary to use an equation which describes the transport of the 

: vu \ 
tracer in terms of the mean (temporal) values u, and c. If the 


flow is incompressible, Equation 2.1 reduces to: 


dc 0 0 dc 
= +-z—(u,c) = -- @ = (2.2) 
ot OX, ft: 4; 2 OX, 


For turbulent flows, we can resolve instantaneous values of u, and 


¢ into mean and fluctuating components such that: 


u, = A +2 ual! 
i i i 
ef (203) 
eased ee eas i 
vv ou _ " d " th 
where u; and c are the mean components; u," and c" are the 


fluctuating components. 


Substituting Equation 2.3 into Equation 2.2, applying the 


Reynolds averaging procedure and using the continuity principle results 


in the following equation: 
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The magnitude of the molecular diffusion term is usually much less than 
the turbulent term on the RHS of Equation 2.4 and can be ignored in the 
present context. The velocity-concentration covariance terms ue), 
represent the transport of tracer by the turbulent fluctuations. Trans- 
formation of Equation 2.4 into a form amenable to solution hinges on 


Wy 


the interpretation of the u, ¢ terms. 


A turbulent diffusion* coefficient, € , is generally introduced 


at this point, just as Boussinesq did for the transport of momentum, such 


that: 
ay 
trictanes ters (2.5a) 
5H OX, 


% Sad ‘ 
Because ute" and dc/ox, are vector quantities, the turbulent dif- 


fusion coefficient, € , must either be a scalar or a tensor of the 
second order (Hinze, 1959). If e€ is considered to be a scalar then 
Equation 2.5a applies. It may be remarked here that it is in general 
unlikely that the turbulent diffusion coefficient will be constant 
throughout the entire flow field. If e¢ is a second order tensor, then 
Equation 2.5a should read: 

Neer e rena SeTEaSUEaT GEERT Ser URES Me ere a Pena ST 
* The terms diffusion, dispersion and mixing as used in this report have 
the following meaning. Diffusion connotes transport that is associated 
primarily with the mean (temporal) product of velocity and concentration 
as defined in Equation 2.6. Dispersion is transport associated with the 


product of mean concentration and velocity differences. Mixing denotes 
transport caused by both diffusion and dispersion. 
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According to Fischer (1970), whether one considers the diffusion coef- 
ficient to be a scalar or tensor is usually of little import, for nei- 
ther is an exact representation of the turbulent diffusion process. 

In what follows, therefore, we shall refer only to the longitudinal, 
Vast and transverse diffusion coefficients, ey? ey Sie = and 
these coefficients will be regarded as local properties of the turbu- 


lence field. 


Apart from the work of McQuivey and Keefer (1972) the velocity- 
concentration covariance terms in Equation 2.4 have not been success- 
cully measured. McQuivey and Keefer measured the longitudinal velocity- 

UYU 


concentration covariance term uc in a shear flow and evaluated the 


longitudinal diffusion coefficient, Ey » from the relation: 


They obtained E. values which were within + 20% of the corresponding 
surface values determined from floating particle studies. They found the 
Boussinesq-type turbulent exchange coefficient to be an adequate model 


for describing the diffusion process in the longitudinal direction. 


Equation 2.4 is now expanded to conventional rectangular carte- 
sian notation with Xp» Xo» ¥3 = Xs Vs 2 and Uj,» Up, Us = Us Vs Wi 1k 


is noted that: 
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Neglecting molecular diffusion and dropping the time-averaging bars we 
obtain the general convective-diffusion equation for a tracer in a 


turbulent flow: 


deme d(uc) (vc) we) Ase oo dc. 
ene dy are Ox — ax? i ; 
(237) 
iy ub. des 5, Sangre 
oy (e, ay? y az ‘“z dz 


In applications to open channel flows, x, y and z will represent 
distances in the longitudinal, vertical and transverse directions res- 
pectively, with u,v and w as the corresponding velocities. 

Equation 2.7 is still too complex for a general solution in 
a realistic situation. Generally, solution of Equation 2.7 for open 
channel flows must account for a variety of factors such as velocity 
gradients in all directions, channel shape, roughness, obstructions, 
bends and the variation in the turbulent diffusion and mixing coef- 
ficients. No existing theory does this adequately. Most investiga- 
tions have been devoted to one or more of the following: 


4s Establishment of some basis for predicting the turbulent 
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mixing coefficients and the concentration distribution 


downstream from a tracer source, 


2. Transformation of Equation 2.7 into forms more amenable 
to solution either by making appropriate assumptions or 


by spatial averaging procedures, 


3. Finding approximate solutions by finite difference 


techniques using digital computers. 


Considerable progress has been made for the simple idealized 
case of uniform 2-dimensional flows which are approximated in labora- 
tory flumes. Investigations of this type of flow have been made by 


Glover (1964), Elder (1959), Sayre and Chang (1968), Okoye (1970). 


2.2 Solutions of the Convective-Diffusion Equation for 
Transverse Diffusion in Straight Channels 


To facilitate solution of Equation 2.7 for uniform turbulent 
flow in a straight and wide rectangular channel, complete analogy with 
molecular diffusion is assumed. The transfer coefficients Es = and 
€, are assumed constant and the longitudinal velocity u is replaced 
by Ue ,» the cross-sectional average velocity. It is assumed that any 
variations of u within the cross-section and any secondary flow 


effects can be absorbed into EL: = and 2 With these assumptions 


Equation 2.7 reduces to: 


dc dc dc r) Cc OuG (2.8) 


Solutions of Equation 2.8 for various initial and boundary conditions 


ie) ee cs 
a9 : 


sotguds t3ebb aol irxan9 snes sid) here 
ie SOTO a a" 


esidansma s7om nie er. ae nas an to 


esaerel tlh stiankt ‘ ine 9s 


gr) i 


bexllasbk ofgake ada 208 Bae saat “ead, sore 
~nrods f nt bodtimbrongds. at prt avolt 


vd obam asd) ove wol} ic aneaigihis ia 
ores) exadd cea)? duet ‘bas ane eine 


— 


tan: ‘yt ty 3 dinetoitigos xohenexa rig ‘bomen, ar fod at 
- leona et v sokpotae sasha tarot ets ‘bate sons mt 
(as “asaa ‘fbmeen ok a \eatsotey eger9ve Fgaprdsbernaaye § ty ah 

wot 3 KaBbAGoNe - ents. hae aobabge~esors of pidgiw a te ana i 
Bnoksqmack vas h2¥ * :? 5 ey on nga) oda, baiiroadh 9d 


I 


’ 
ay 


v 


02 ssouhen NS aokanopl 


j : 4 im mys’ 
(8.8) | a cies age, 2 ed. oueh' wn 


eutadindinn ‘inbred bats inti dione 163-8. eorsauel io suotsubet 


ae 


have been presented by Carslaw and Jaeger (1959), Glover (1964). For 

a continuous point source in a straight peeraneulat channel with a 

large width-depth ratio, Equation 2.8 can be reduced to 2-dimensions 

by integrating with respect to depth and noting that there is no tracer 

transport across the wetted perimeter and across the water surface. 

Equation 2.8 then becomes: 

© a erhce ae oc (2.9) 
Ox 

where c = c(z,x) is the depth-averaged concentration. The solution 


is (Sayre and Chang, 1968): 


Ux 
fe) 
ae Q_c ae U or Ries 
c(z,x) = oe tere te tpl scion, Safe (2.10) 
2nd o}2e S 
® AE oe x Zz 
in which Q, = volumetric tracer discharge 

C. = tracer concentration at source 
d, = flow depth 


modified Bessel function of the second kind, 


A 
i} 


of order zero 


When the distance downstream from the tracer source is large enough 


such that: 


€ 2 2€.. 
(2)5 <<a of and x Sp == 


x Oo 


Equation 2.10 converges to: 
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c(z,x) = —_—__——— (2.11) 


This represents the solution of Equation 2.9 under conditions of neg- 
ligible longitudinal diffusion. Equation 2.11 can be applied to a 
straight channel of finite width provided the tracer is well mixed 
over the depth and a significant amount of tracer has not reached the 


side-walls. 


A more general solution of Equation 2.9 applicable to channels 
of both infinite and finite width is given by ee and Chang (1968). 
For the general case where the tracer source is located a lateral dis- 
tance from the origin (center of channel), they employed the method of 


superposition and treated the side-walls as reflecting barriers to 


obtain: 
co 
c(z,x) = c(z-T,x) + © [ctnW- TF + (-1)"2,x) + 
n=1 
se A (2512) 
c(- nW - t + (-1) z,x)] 
where . -W/2 = ia, = W/2 and n = number of reflection cycles. The 


terms on the right hand side of Equation 2.12 are evaluated from Equa- 
tion 2.11. Equations 2.11 and 2.12 are in fairly good agreement with 


experimental observations in straight rectangular channels. 


Another method for solving Equation 2.7 was recently proposed 
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by Yotsukura and Cobb (1972) for straight prismatic channels of arbi- 
trary cross-section. They neglected longitudinal turbulent diffusion, 
assumed a steady state, but accounted for transverse variation in depth. 
They showed that the transverse diffusion equation for this situation 


is: 
wet = 2 Ens (2.13) 
where u, D, Ce are functions of z only; eo and u_ being depth- 


averaged values. If the origin is located at one bank, the boundary 


condition at the channel edges (z = 0, W) is: 


—-. de _ 
ED ope Ow 0 (2.14) 


The conservation of tracer is expressed as: 


q 
M 


W 
y[aevaz (2515) 
ie} 


where y is the specific weight of water. 


Yotsukura and Cobb then defined a new independent variable 


q,» the cumulative discharge, measured from one bank as: 
Zz 
do. ox | uDdz 
fo) 


This can be used as an alternative transverse co-ordinate. By using 


q, to replace z they transformed Equations 2.15, eels and 2.15° to: 
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dc 9 — —2 9¢ 
qe fer Sggeesad san page: (2.16) 
- —29c _ 
S D 3a, = 0 at q,-= 0) -andy iQ 
and: . aye = ¥ cdq, 
fe) 


respectively. To obtain a solution to Equation 2.16 they assumed that 


it can be approximated by: 


es oe 
per Locale — 2omdee 
aot E uD 9 (2-16a) 
oq, 
: . ee Wee 
where the Diffusion Factor € uD = Q E uD dq,- 
° 


They finally applied the method of images used by Sayre and Chang’ (1968) 
to obtain a series solution for Equation 2.16a. This theoretical model 


seems to agree well with field data. The approximation of the variable 


qe ; ; 
parameter € uD by its cross-sectional average E,ud however needs 


more verification. 


2.3 Evaluation of the Transverse Diffusion Coefficient 


for Straight Channels 


The turbulence which is responsible for vertical diffusion 


also causes diffusion across the channel. The vertical diffusion co- 


| 


(a2.2). 
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efficient can be derived by assuming (i) equivalence between mass 
and momentum transfer (Reynolds analogy) and (ii) that the shear 
stress vee =- ony is linearly distributed over the flow depth. 


The following relation is obtained: 


uy 


AV) 
MYUY ee (tie 29/0) 
€ in Settee Waments (2.17) 
7 du du 
dy dy 


where T is the bed shear stress. So far, no analytical relationship 
has been derived between the vertical and transverse diffusion coef- 
ficients. The transverse diffusion coefficient, ey » Must at present 
be evaluated empirically from experimentally observed concentration 


distributions. 


The local transverse diffusion coefficient, ee » in a two- 
dimensional open channel flow is known to vary over the flow depth, 
(Okoye 1970, Fischer 1973). It increases from zero on the channel bed 
to a maximum near the free water surface. In general, one is interes- 
ted in the overall spreading of a contaminant in a stream and the use 
of an average transverse diffusion coefficient over a channel reach, 


E. » simplifies the computations considerably. 


Two approaches have been used in the past to obtain E. from 
the convective-diffusion equation. The first approach is the Method 
of Moments first introduced by Aris (1956). Basically, the method in- 
Pp 


volves multiplying the governing convective-diffusion equation by z 


and integrating with respect to z from -©® to ®. The resulting 
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equation is then expressed in terms of the moments of the experimental 


or theoretical c distribution which are defined as: 


co 


m, (x) = | zP’c(z,x)dz 3; p=0, 1, 2... (2.18) 


—- © 


and then solved for the various moments. 


-In principle, any desired degree of detail can be achieved by 
solving for high values of p. In practice, however, the zeroth, 
first, second and third moments are sufficient for calculating statis- 
tical parameters which are useful for describing some aspects of the 
diffusion process. The method of moments has been used extensively, 
for example, by Okoye (1970), Sayre and Chang (1968), Holley et al 


(1972). 


For a continuous point source in a 2-dimensional open channel 


flow, it can be shown by the method of moments that: 


ue do,” 
a ea0. 2.19 
= ee ony 


provided the tracer cloud has not encountered the channel walls and the 
locus of the centroid of the transverse tracer concentration distribu- 
tion at various sections downstream of the source is parallel to the 
channel axis (z=0). In Equation 2.19, ae = m,/m and is often 
called the variance of the re distribution; Mm, and mM, > for this 


situation, are the zeroth and second central moments respectively. 
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Equation 2.19 can also be obtained by first finding the vari- 


ances corresponding to Equation 2.10 and Equation 2.11 which are: 


2 2E. 2E. 
Toss geet 7G (2.20) 
(0) ° 
2 2E x 
and PR u, (221) 


respectively. Differentiation of Equations 2.20 and 2.21 then gives 
the same result as Fauation 2.19. The second method used for finding 
ED is designated "the simulation method" (Yotsukura et al, 1970). It 
is a trial and error substitution of the appropriate diffusion coef- 
ficient into a numerical or analytical solution of the governing equa- 
tion until the best fit between calculated and observed concentration 


distributions is achieved. 


A summary of previous measurements of the transverse diffu- 
sion coefficient in open channel flows is given in Table 2.1. Labora- 
tory experiments with dissolved tracers in straight flumes have yielded 
values of the normalized diffusion coefficient E/u,d,, (=k) ranging 
from 0.10 to 0.23, where u, is the shear velocity, d, is the 
overall average depth, and E.. is the average value of transverse 
diffusion coefficient over a channel reach. Okoye (1970) found the 
aspect ratio, defined as W/d, ,» to be the most significant parameter 
controlling the transverse spreading of a tracer injected continuously 


into a straight rectangular open channel flow. The dimensionless 
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diffusion coefficient, k , increases with an increase in the aspect 


ratio. 


2.4 Transverse Mixing in Meandering Channels 


The characteristics of transverse mixing in turbulent open 


channel flows with curved boundaries is extremely complex. The Fickian- 


type Equation 2.9 which successfully predicted the spreading of tracer 
eee straight uniform channels with a known E, is found to give un- 
satisfactory results for curved channels. This is attributable to a 
combination of various factors such as depth, width and velocity 
variations, growth, decay and reversal of the secondary motion and 
variable turbulence structure. Recent studies by Fischer (1969), 
Yotsukura et al (1970) and Chang (1971) have, however, contributed im- 


mensely to our understanding of the mixing in meandering channels. 


The first study of transverse mixing in bends of turbulent 
open channel flows was reported by Fischer (1969). He derived an ex- 
pression for predicting the lateral dispersion coefficient, Ea» based 
on the following assumptions: (i) lateral dispersion is analogous to 
longitudinal dispersion and (ii) that the transverse velocity distri- 
Piciee ti the vertical can be described by Rozovskii's transverse 
velocity distribution for fully developed turbulent flow in an open 
channel bend. Fischer verified his analytical findings experimentally 
in a meandering flume having a 330° bend with a radius of curvature of 


206 cm and a width of 76.3 cm. His observed values of Ea agreed with 
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the predicted values within a factor of 2. The large discrepancy arises 
because the bend was too short for the secondary flow pattern to have 
reached the fully developed stage. Fischer also attributes the errors 
to the sensitivity of the result to small changes in the input para- 
meters, which are often not accurately known. The procedure is not a 
practical method for estimating the transverse dispersion coefficient 
in natural streams because the condition of fully developed bend flow 
is never obtained in natural streams. Nevertheless, Fischer seneers 
ted that the dominant effect of curvature is generally to increase the 
normalized lateral mixing coefficient compared to that usually observed 
in a straight channel. His normalized dispersion coefficients ranged 
from 0.52 to 2.4, and this is about 2 to 10 times larger than the ob- 


served values for straight laboratory channels. 


A simulation procedure was used by Yotsukura, Fischer and 
Sayre (1970) to predict the transverse mixing coefficient in a six-mile 
reach of the Missouri River near Blair, Nebraska. The procedure is a 
numerical solution of the convective-diffusion equation. The steps 
involved are briefly explained. The river channel is first divided 
transversely into 20 stream tubes carrying constant and equal discharge. 
The mass transfer between stream tubes is assumed to be entirely diffu- 
sive and any effects of secondary currents are absorbed in the diffu- 
sion term. The governing equations are then expressed in finite 
difference form and solved by trial and error until the mixing coeffi- 


cient gives a calculated concentration distribution close to the observed 
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data. For the six-mile reach of the Missouri River near Blair, Yotsukura 
et al obtained an average value for FE. of 0.6 u,d, by the simulation 


technique compared with 0.7 u,d, from the Method of Moments. 


Chang (1971) proposed two methods for finding the lateral 
mixing coefficients, namely, a simulation procedure and an integral 
method. His analytical derivations are based on the convective diffu- 
sion equation expressed in a meandering co-ordinate system (Sayre fas 


Fukuoka, 1973): 


ctimd £9 Cen) a Cev) Bebe 1 te ALL. OO dc. 
se ss h ox as Oy : h Oz 200x (e,. Dx) ui 
1 al hy 

(2.23) 

3 dc 1 39 dc 

Oy ( y ay) i h, dz (hye, Oz 

and the continuity equation: 
0 (hw) 

du ov ie 

= x +———-_ = ae 

a + hy By + De 0 (2.24) 


where the x-axis coincides with the centerline of the channel; y is 
measured downward from the free water surface orthogonally to the x- 
axis, and z is measured laterally from the channel axis orthogonal 
to both x and y. Depending on whether the channel is straight or 
curved, the surface z = constant can be either a plane parallel to 


the x-axis or a cylindrical surface concentric with the x-axis. In 


Equations 2.23 and 2.24: 
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hy = ji + a for a bend curving to the left 

c 
hy = 1 for a straight reach A are) 
hy = l- a for a bend curving to the right 

c 


and te is the radius of curvature of the channel axis. 


The simulation technique employed by Chang is similar to that 
used by Yotsukura et al. He integrated the steady state form of Equa- 
tion 2.23 and Equation 2.24 over the cross-sectional area of the j'th 


stream tube to obtain: 


Zz 
oCe A ae j A 
aoe — = h,D(e, + €,) ae (2.26) 
z.—-1 
j 
in which 2 Be discharge in a stream tube, Z denotes the boundary 


between the (j+l)'th and j'th stream tubes, C, is average concentra- 


tion with the j'th stream tube, Ey is the transverse dispersion coef- 


ficient given by: 


US LOR ae (2.27) 


where w' and c' are deviations from the depth-averaged values of 


w and ec, 


Chang then solved Equation 2.26 numerically to find the best 


- value of the lateral mixing coefficient, E, (= &. + E4) » which gave 


good agreement between the calculated and observed concentration distri- 
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butions. The use of Equation 2.27 is subject to debate since negative 
values of €, can be obtained. In such a case the dispersion ceases 


to be of the gradient type. 


The Integral Method rests on the same reasoning as the pre- 
vious method but provides a direct method for finding the local lateral 
mixing coefficient from experimental data. Chang integrated Equation 
2,23 and Equation 2.24 over an area extending from the left bank to any 


desired stream tube boundary Z4 and arrived at the relation:. 


eieD 
) 
| | uc dy dz 
3 e) Oo 
ee a a |. (2.28) 
by? 5 
ot 


Chang's results indicated a strong periodic variation of the 
lateral mixing coefficient in the longitudinal direction with the 
simulation and the integral methods giving similar patterns in the 
variation of Ew’ However, the integral method predicted values of 
E ow approximately twice as large as the simulation method. It is 
noted that the integral method is more susceptible to errors because 
of the inaccuracies involved in taking derivatives of experimental data 


Chang obtained normalized mixing coefficients ranging from 0.6 to d.2 


with the simulation method. 


Holley, Siemons and Abrahams (1972) have put forward a method 
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for predicting the lateral mixing coefficient, Ew (the Generalized 
Change of Moment Method). Moments of the tracer flux (D uc) distribu- 
tion are used rather than that of the concentration distribution. The 
method incorporates transverse variations in depth, turbulent exchange 
coefficients, longitudinal and transverse velocities into the depth- 
averaged convective diffusion equation. They tried various expressions 
for the variation of the exchange coefficient across the channel. 


These were: 


EO, 7 kdyu, = k,4,U, (2.29) 

Eve) ROS (fa (2:30) 
U 

Uc os Dees (2.31) 

ZW va, 


where d, is the overall average depth and ki ks k, are dimen- 
sionless mixing coefficients and ee is the local depth-averaged 
mixing coefficient. US is used in Equation 2.29 simply as a matter 
of convenience. In Equation 2.30, ce is assumed to vary as the 
product of local depth and local depth-averaged velocity while in Equa- 
tion 2.31 Gus is taken to be proportional to the product of local 


depth and local shear velocity. Equation 2.31 is arrived at by use of 


Chezy's equation. 


Holley et al recognized the limitations of Equations 2729, 


2.30 and 2.31. They found that the type of variation assumed for er 


can have a significant effect on the value of k, or k, determined 
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from measured concentration distributions. Nevertheless, the use of 
local values of depth, velocity and shear velocity appears to be rea- 
sonable since they represent the accumulated effects of various up- 


stream factors which also influence the local value of e 


This chapter has presented a review of previous studies on 
transverse diffusion and mixing in straight and curved open channel 
flows. No laboratory or field measurements of the transverse diffusion 
or mixing characteristics in ice-covered channels appear to have heen 
reported so far. Investigations on transverse diffusion and mixing in 


ice-covered channel flows are therefore urgently needed. 
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CHAPTER 3 


ANALYTICAL INVESTIGATIONS 


3.1 Similarity Analysis for Straight Channels 


Consider the case of a continuous release of tracer at a 
steady rate into a steady uniform turbulent flow in a wide straight 
rectangular channel with no side-wall effects. If longitudinal dif- 
fusion is assumed to be negligible and E. constant within test reach 


then Equation 2.9 reduces to: 


For practical application at can be replaced by E. » the average 


over a reach of length x , such that: 


Aa Ee (3.1) 


To make Equation 3.1 dimensionless, the following dimensionless vari- 


ables are introduced into Equation 3.1: 


U — 
Sh nay aa Co reds 
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where A is the aspect ratio and Ge is tracer concentration when 


25 


25 


6 38 tecex? lo aaadley eoGuRRInAD. a toaens 9 


} ee 


tigiezie abiv x salle Sa aa seh ’ - ne se ; 


-tth fenatbysiehnel 3 attatte ious es gta Lan 
doses geo abdtiw sesteceso * or ‘haiiiee 
“et oni 


a x ¢ cad De: 
ae 3. =e Preion: | ¥ ites 
5 c Bee Tee, oO” ; es a 7 : 

oe ‘se fae . ine 


egeteves sft. 3 yd beasiqed 'ad mg"! < The: 


46rd dove, % dagast Yo is 


| 


“aay sasinotensiitt sitvotied its _ seated ienbit £8 poner ja 
vine we naan tnt: Scouboraet oa a 


. 
oe 
: Raly 


: ‘ . + ee 


i Weal ind £m SY 
) a pete: - an co 
. ae " i 5 
ie , Ale re an ellen a Ue We he 


26 


fully diluted. The resulting equation is: 

ac arent 
= k— (3.2) 
1 es) 
It is a well-known fact that as the tracer cloud oe the depth- 
averaged concentration along the tracer centerline decreases. Experi- 
mental evidence (to be given in Chapter 6) indicates that the normalized 
concentration profiles, Cy (x,8) , at various sections can be made 
congruent by appropriate choice of scaling parameters for Cy and 8. 
Under such conditions the profiles are said to be "similar". The re- 
quirement for "similar" Cy profiles may be stated as the property that 
any two C) profiles located at different X differ only by a scale 


Factor in C and §. This can be stated mathematically as: 


2 


B B 
c (X_) ——— Cc (X, ) Da) ipa areas 
Ly jet ai 1 2 gl (X,)5] 
o Ain and Pe ls Nyg®  volieeian 48 new 2 vou ti(ong) 
Cin [(X,),] Chan [(X,)5] 


where Cyn OX 


respectively and Ck is the normalized peak concentration. 


) and g (X,) are the scaling parameters for C, and 8 


An obvious characteristic scaling parameter for the tracer 
cloud would be its dimensionless width at the section under considera- 
tion. However, exact demarcation of the cloud boundary would be dif- 
ficult. An alternative scaling parameter is the dimensionless trans- 
verse distance between the tracer centerline and the vertical along 


b 
= = r B 
which C, (X)> B) 4 Fo) Cy OX , denoted here by b, &) ( W ) 
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Note that b is one-half the transverse distance between the verticals 


along which C=0.5 Ci and Ci is depth-averaged peak concentration. 


If similarity does exist, one could assume the following re- 


lation: 
al 
ae et (by b-) = vt (a) (3.4) 
Gc 1 
im 
where Qa = B/b,- Substitution of Equation 3.4 into Equation 3.2 
yields: 
2 
e) | er) 
ea Ce.f =) eam a TCE] 
dX, im U 982 1m 
f or Coe! ea ay ae “1m ¢" (3.5) 
OX, 1m by dX, Uy oe 


2 


where f£' = df/da and f" = d7£/do . Equation 3.5 can be rewritten 


as: 


ce ae keeetim aeOf. me aby k_ “1m Ee hy (3.6) 
£1 3X U Z oy babe § ak U, b 2 ‘ 


The R.H.S. of Equation 3.5 is independent of x » therefore the L.H.S. 


must also be independent of Xj This implies that: 


oC G 
Fale “s| = constant « ue (3.7) 
ib i by 


c db: G 
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If it is assumed further that: 


and: (3.9) 


where p and q are as yet undetermined exponents, Equation 3.6 or 


Equation 3.7 yields: 


(3.10) 


a 
" 
NS) le 


To evaluate the other exponent, p , an additional equation is required. 
This is obtained from the condition that there is no loss or decay of 


tracer, that is: 


on 4,u, | edz (3.11) 


where is » the rate of injection of tracer into the flow, is equal to 
UA: C.. A, is the area of the tracer injector and S is concentra- 


oOo Ss 


tion of tracer at source. 


When Equation 3.11 is simplified using the appropriate scaling 


parameters, the following relation is obtained: 


Ag(C./C.) 
b.C Sy Sed (3712) 


2 
ae | fda 


The R.H.S. of Equation 3.12 is again independent of x, » so that the 


L.H.S. must also be independent of x) for the equation to be valid at 
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X q xX P = constant © Xy 


ap ee. 0 (3515) 


L 
2) =e (3.14) 
Therefore: 

A 

Caley 1/2 oe Cm ub (3.15) 
VE) 

Ly2 ee 
bo, or bs AVX, (3.16) 


where Ay and A, are constants for a particular flow situation. The 
form of the function f may be found by first substituting Equation 


3.15 and Equation 3.16 into Equation 3.5. The following equation is 


obtained: 
a £"OOSACGE*. +ANE Ci atolO 
UA, 
or: 
yim fF" 4+ © (af) = 0 Gals) 
UA da 
i ay 


Equation 3.17 can be integrated twice with respect to a , recognizing 


that when a = 0, f(0) = 1 and f'(o) = 0. The solution of Equation 
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3.17 becomes: 


f(a) = e = — (3.18) 


The distribution of C,/C, is therefore Gaussian. Equation 
3.18 also implies that for a given flow situation in a very wide rec- 
tangular channel, away from side-wall effects, the concentration pro- 


files at various downstream sections is described by a single curve. 


When a@=1,C, = 0.5 C = and from Equation 3.18 we get: 


ba ee eee (3.19) 


Combining Equation 3.19 with Equation 3.16 we obtain an expression for 
E_: 


bu 
e) 


z 4xi1n2 (3-20) 


Equation 3.19 suggests that the normalized transverse exchange or dif- 
fusion coefficient depends only on UL (= Uo/Y») or the friction factor 
ne A, is a universal constant. The effect of width was effectively 
removed by confining the analysis to the central portion of the flow 
away from the side-walls. Although most of the above results are well- 
known, the derivation is novel and shows that for a given flow situation 
in a rectangular channel the lateral concentration profiles at various 


sections can be represented by a single curve. Equation 3.20 also pro- 


vides a quick method for evaluating the transverse exchange coefficient, 
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E. » given the average velocity Ue and the lateral depth-averaged 


concentration profile at a downstream section. 
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The Two-Dimensional Lateral Diffusion Equation 


The convective diffusion equation for turbulent flow in a 


meandering co-ordinate system may be written as: 


proc 8c) 0 (vc) Os ~ 
Ge obyeiox, jondveeMbani) 22fsn2ice oo 
(2523) 
FINS fe) uea ek gc ae 
h 2 ox (e% 3x) i oy So ay? of hy Oz (hye, oz 
1 
where - hy, = 1+ 2/r_ for a bend curving to the left 
= 1l-2/r_ for a bend curving to the right 
= 1 for a straight reach 
EL» ey? as are the turbulent diffusion coefficients in the x 


tudinal), y (vertical) and z 
the time-averaged velocity components in the 


respectively; and c 


(lateral) directions; 


x50 = and. 2 


is the local time-averaged concentration. 


(longi- 
U, Vv, W are 


directions 


There is no mass transfer across the channel boundaries or 


the water surface by either convection or diffusion. 


conditions can be stated mathematically as: 


8 ioe a oc QO at the water surface 
y oy 


w= Ossev =03 
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boundaries 


These boundary 


(3221) 


= 0 at the solid 
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Solution of Equation 2.23 requires a detailed knowledge of the distri- 
butions of local velocity and turbulent diffusion coefficients and it 
is also very cumbersome to work with Equation 2.23. To facilitate 

the analysis, the equations can be reduced to 2-dimensions by averaging 


over the local depth, D. 


Consider non-uniform flow in a channel of arbitrary cross- 
section and a continuous, steady release of tracer into the channel at 
a fixed location. Longitudinal diffusion does not have a significant 
effect on the steady-state concentration distribution and can be neg- 
lected (Sayre and Chang, 1968). If the transverse diffusion coeffi- 
cient, Ee. » 1s assumed not to vary AA AS over a considerable 
portion of the flow depth then the steady-state case of Equation 2.23 


can be written as: 


1 (uc) d(wve) Ld 2 
— +—-7 (h,wce) = 
hy ox dy hy vere 1 
(8% 22) 
eee eee: 
ay ay i hy dz caret 32> 
where os is depth-averaged transverse diffusion coefficient. 
Define: 
u(x, y, Zz) = u(x, z) + u'(x, ys 2 
yG@ityiaz) c= v(x, z) + v'(x, y, 2) 
; (3.23) 
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where the bar again denotes a depth-averaged quantity (for example, 
D 
u == Jf udy) and the prime designates a local deviation from the 


2 ° 
depth-averaged value. Equations 3-23 are now substituted into 3.22 and 
the resulting equation integrated term by term with respect to y from 
zero to OD. The order of integration and differentiation are inter- 


changed, where necessary, using Leibnitz's Rule. The resulting 


equation is: 


Opr@rt caste te—F: on. BOs ee a ee 
ox we + uc yO a {Dh, (w c+w'ec')} = 
Oy (3.24) 
P] — oc 
az tPhYE, oz 


The first term in Equation 3.24 refers to the longitudinal differential 
convection or dispersion. Except close to the tracer source the con- 
vection dispersion term utc’ is usually much smaller than uc. and 
can often be neglected. The second term represents differential con- 


vective mass transport in the lateral direction. If a Us ow, cu 


tA! 


Cc 
and w'c' are known or can be predicted, Equation 3.24 can be solved 
for the steady-state concentration distribution at any downstream sec- 


tion from the tracer source. 


Normally, eo is assumed to have a constant value within a 
channel reach but it seems probable that significant changes in depth 
or velocity within a channel reach may cause significant changes in é.- 
At present there is, however, insufficient information about the turbu- 


lence responsible for diffusion to permit mathematical representation 
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of the variation of €. For the time being, the variation of ce 


can be expressed as follows: 


ia = E!o(x, z) (3.25) 


where E. is a constant for a given flow condition and @$ defines the 
variation of Cm It is to be noted that when ¢=1, E = E. » the 
average transverse diffusion coefficient over a channel reach. As 
mentioned previously in Section 2.4, some plausible assumptions regard- 
ing $ have been put forward by Holley et al (1972). However, no 
matter what assumption is made for 96 , E. has to be evaluated ex- 


perimentally. Introducing Equation 3.25 into Equation 3.24 gives the 


2-dimensional lateral diffusion equation: 


ie {D(uc tu'c')} + {Dh, (wc +w'e')} = 
= (3. 24a) 
-) dc 
EY ae {Dh,¢ aa 


Since E, is the principal transverse diffusion parameter subsequent 


analysis and discussion of results particularly related to the straight 


flume experiments will concentrate on it. 


3.3 Methods for Evaluating the Transverse Exchange Coefficient 


Three general methods are presently used to evaluate the 
transverse exchange (diffusion or mixing) coefficient, E. » in channel 


flows. These are (i) The Simulation Method (ii) The Method of 


34 


{3% ootestisy om cgister ens ae tot 


‘, 


= dyes Re tis 


Ae re a 


(25.0) | ere a) 


oi2 asnkish 4 Sas cotsibimo wolt aoety 8 fo 4 a | et 
2. f= . fs 6 ite onda wane ar 
aA 6.isees Lontigde 9 xaVve sneksthises oka 3: Ke 
-ipget eoriqmgens. sftyeusle, OFA BS) actin ghia 
on ,wevewol . CONE). dp 39 vation xd hisiwro2 2004 
-«4 bespuleve ad ot) eat: 7 ane 102 obi" et b3 
ads aevig 6S.£ colseupal otnt Wt aOF 2013, & sub 


(ab. ) 


iG w = a 
sasupoedue iw’ ‘ona8@ naeutite sexoreners teatoatig “4s “a of — ; 
a 


shabarsive aaa Jos bedalsx Mratvosacbq etirees. to sstonvontt ‘Rea carci | 
' + 9 ‘Aes 
at ae Faia ‘tt seat oayth 


‘parae ay 


ad? eagutave of bean isnaagig bev wi i sett Noe 
| | Ai pled yey, 58 
" Eooneds ah. 2 .sngkokiedy Contin 20 haw) piiox 

Bo: hottett ott ce beabit aoksitun sat a se west? seve 


35 


Moments and (iii) The Integral Method. These methods were briefly 
reviewed in Chapter 2. Essentially the same amount of tracer concentra- 
tion and hydraulic data is required by all the three methods for eva- 
luation of E. » but the Integral and Moments methods have the advantage 
of computational simplicity. On the other hand, the procedure of 
simulating the mixing process is superior to the Integral and Moments 
techniques because it can be used to predict the tracer concentration 
distribution downstream of a tracer source if the transverse exchange 


coefficient, E. é is known. 


Two methods for evaluating E. based on the Integral and 
Moments methods are proposed here. The method of moments is used in 
the present analysis to determine an overall transverse exchange co- 
efficient within a channel reach. The Integral approach attempts to 
determine the distribution of the average transverse exchange coef- 


ficient across and along the channel. 


3.3.1 Method of Moments 


The method of moments was first described by Aris (1956) and 
has since been used widely in numerical work. It suffers from the 
disadvantage that no direct expression for either the ce or Duc 
distribution is obtainable. However, knowledge of the first few moments 
of the c or Duc distribution gives a great deal of information 
about the distribution itself. Although in principle higher order 


moments of a distribution about any value can be computed, a meaning 
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cannot necessarily be attached to the resulting values. In practice, 
only the first four moments are used because accuracy in estimating 
higher order moments from experimental data decreases rapidly as the 


order of the moments increases. 


In applying the moments technique here, Equation 3.24a is 


2 
multiplied through by n and then integrated with respect to Nn to 


give: 
Wo 
) { t t } 2 
Ox Dtuc + u'c') n dn + 
ui 
W 
; Pe) Wes Sl es ae Zz 
| creme! = Y(wetw'e')} nodn = 
W . , 
1 
W 
2 
ipl Sp — 
P) o oc 2 
‘ eS +—— TE ~26 
Et | mn {D(1 r ) o an? n dn (3.26) 
us) 
Ee Wi) ar v4 
in which N=2z- Zz, 3 hy =l+t ro bee 5 es zy represents 


the transverse distance from the channel axis to a fixed but arbitrary 
cylindral surface which is concentric with the channel axis. For the 
present analysis Z5 is taken as the lateral distance from the channel 
axis to the cylindral surface passing through the point of tracer 
release and concentric with the channel axis. Note that in Equation 


3.20; Wy and Wo represent transverse distances to the channel edges 
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The second and third terms in Equation 3.26 can be integrated 
by parts. Some of the resulting terms go to zero because there is no 
transfer of mass at the flow boundaries. Again, Liebnitz's rule for 
interchanging the order of integration and differentiation is used. 
With some further rearrangement of the z-diffusion term and after 


dividing through by the tracer flux which is constant, we obtain: 


W W 


RN ee : 2 teen eens 
| D(u c +t u'c') n dn 2 | Dn (1 + pores) \ ec + w'e')dn 
(e 
ele ek pes 
dx 9 Wo 
foucan foucan 
Wy Wy 
W W 
2 2n + Zo F he not Zo 3D 
| wea r yan+ | gen ae 2 ) oy an 
W 
iy Sects 1 
2 
[ pucan 
7 
W 
; 3 — Need és 
| » an {Den (1 + ee dn 
W 
~ 27 | - (3.27) 
Zz 
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The diffusion term is found to be a combination of three terms. Term 
A is expected to remain reasonably constant along the channel. Term 
B expresses the Eubertanee of depth variations within a cross-section 
and goes to zero when D has a constant value across the section. 

Term C contributes to the diffusion term only when a significant 
amount of tracer reaches the side-walls. The magnitude and longitudi- 
nal variation of the different terms in Equation 3.27 will be discussed 


in Chapter 6, 


Equation 3.27 is proposed and used in this report for evalua- 


tion of ES from known distributions of w, u, c, u'c', w'e', D and 
d. The above analysis differs from that of Holley et al (1972) only 
in the use of the meandering co-ordinate system but illustrates more 


clearly the effects of curvature and depth variations on the mixing. 


One limitation in using simplified forms of Equation 3.27 is 
worth mentioning. In principle, the same value for E will be ob- 


tained irrespective of where mn is measured from provided the 


variations of u, w, c, u'c', w'c', D and @$ within the reach under 


consideration are correctly accounted for. If 6 does vary within a 


channel reach, as is generally the case, but is assumed constant (¢ = 1), 


then Equation 3.27 will give different values of E depending on where 


1) is measured from. In such a case E. (= E,) can be regarded only 


as a first estimate. Discussion of this problem is deferred until 


Chapter 7. Equation 3.27 can now be applied to several flow situations. 
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Case 1 - Meandering Channels 


If the distribution of w is unknown, which is generally the 


case, the effects of transverse velocities on the mixing must be ab- 


sorbed into a modified exchange coefficient, el. From Equation 3.24 


we have: 
nn Gate ame ae oo sn, (&. oS ea aula) 
ox Oz ieee OZ 
€ is defined ase ac =€ Be -we-w'c! In analogy with the 
zw Zwoz Zz OZ BY 


earlier definition of E ‘ om can be written as: 


€ = = 4 G, 2) (3.28) 


where Bg is constant within a channel reach for a particular flow 


Equation 3.27 can 


condition and = describes the variation of = 


now be simplified to give: 
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Equation 3.29 is well suited for application to both laboratory and 


field data. It is to be noted that since D = 0 when fn “WwW and 


for most rivers, the last term in the numerator of the RHS of 


= 


Equation 3.29 goes to zero. Equation 3.29 then assumes the form: 


Wo 


pads Steryncan 


al la a 
dx Wo 
[oucan 
Wa 
W W 
[ope ee at Joana er es 
Cc Cc 
W W 
2E' 1 
zw Wo (3230) 
| Ducd n 
We 


Equation 3.30 provides a direct means of evaluating the transverse ex- 


change coefficient, ie , if the lateral distributions of u, c, oe 
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and D within a number of cross-sections are known and a reasonable 
assumption can be made for do Subsequent discussions of results 
related to the meandering flume and field investigations will concen- 


trate on “E" . 
zw 


In a practical situation, the meandering co-ordinates may be 
fitted to a natural river using a large scale map of the river. The 
thalweg is taken as the longitudinal, x , axis. If the thalweg is not 
(PREBE known, the centerline of the river at the water surface or 
ice-cover underside may also be taken as the longitudinal axis. The 
river reach of interest is divided by trial and error into sub-reaches 
which are either single bends that can be approximated with circular 
arcs or straight reaches. The centerline radius of the bends are then 


geometrically determined. 


Transverse distances are measured normal from the x-axis in 
the horizontal plane. these distances are considered negative when 
measured towards the center of curvature and vice versa. Vertical 


distances are measured downwards from the free water surface or ice- 


cover underside (see Figure 3.1). 


Case 2 - Straight Channels 


For straight channels with arbitrary cross-section tae oo 
but w is not necessarily zero. Transverse velocities caused by 
changes in channel geometry with longitudinal distance may significantly 


alter the mixing. If transverse velocities can be accurately determined 
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then the relevant equation to use is Equation 3.27 with swans Other- 
wise the effects of transverse velocities must again be included in the 
diffusion coefficient as was done in Case 1. The appropriate equation 


in this case is either Equation 3.29 or Equation 3.30 with eee 00, 


For the simplest case of uniform turbulent flow in a prisma- 
tic rectangular channel 9D/3n = 0. If u'c' is negligible, Equation 


3.27 reduces to: 


Wo Wo 
cpsaariatiie 2E! pret. Shy 
z on 
aver Wy 
dx | W = W (3.31) 
Se 2 
fucan fucan 
Wa a 


Away from the side-wall boundary layers u is a constant (= UY and 
@ can be taken as equal to 1.0. If Ed is replaced by E, and 
provided no tracer has reached the channel walls, Equation 3.31 can be 


written as: 
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Equation 3.32 is the well-known formula for finding the diffusion co- 
efficient for uniform turbulent flow in straight rectangular channels 
(Sayre and Chang, 1968, Okoye, 1970). Going through the various flow 
situations has been instructive in demonstrating the general applica- 


bility of Equation 3.27. 


3.3.2 Integral Method 


Integration of Equation 3.24a with respect to z , from one 


bank (denoted by W,) to an arbitrary z _, gives: 


z Zz 
[& {D(u c + u'c')} dz + [3 {Dh, (wc + w'e')} dz = 


Wy Wy 


Zz 
3 dc 
v a eels 
ae | OZ {Dh,¢ 32) ae 
a 


which can be simplified to: 


zZ — 
Sk hese mtn wee. 5 
a | D(u c teu'c')dz + Dh, (w ec +tw'e') = E Dh, oe (3.33) 
Wy im 
=- c 
SEs Dh, ae (3.33a) 


Equation 3.33a can be used to determine the lateral variation of €. 
across a channel reach. Chang (1971) used a similar procedure but di- 


vided the channel into stream tubes carrying equal discharge. His 
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z,-values represented the boundary between the stream tubes and he con- 


bined the dispersion term w'c' with ee dc/dz by assuming w'c' = 
Ey dc/dz where Ey is a dispersion coefficient. The only drawback 
in using Equation 3.33a is that it is susceptible to errors peeatiaetor 
the inaccuracies involved in taking derivatives of experimental tracer 
concentration data. To circumvent this difficulty, Equation 3.33a can 


be integrated once more with respect to z from My to z , to give: 


z Z z 
a> cae ae Br ET as =e 
| 1% | D(uctutc daz} dz + | Dh, (w ec + w'c')dz 
al Li By 
Z — 
oc 
~ ' ee: 3 
= E} dl Dh,¢ a dz (3.34) 
ak 


Without too much loss in detail an average $ can be assumed between 


= E! b ation can be 
Wy and z_ such that E. eve and the above equ 
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Equation 3.35 will give some indication of the variation of the aver- 


age exchange coefficient over transverse portions of a channel reach by 
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performing the above integrations for different z-values across the 
channel. If the effects of transverse velocities are not separated 


but are lumped with the diffusion term in the usual way, then Equation 


3.35 can be modified to give: 


z Zz 
de coe sa 
| ag | D(uctu'tc aed 
bal oY 
Ura 7 — 9@h,) 
= E ow Dhjc - | c — 7, a2 (3.36) 
Wy Wy 


Equation 3.36 is a new expression which is proposed for evaluating ES 
provided variations of u, c, u'c' and D within a channel 


reach are known or can be estimated. 
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CHAPTER 4 


LABORATORY EQUIPMENT AND METHODS 


4.1 Introduction and Experimental Objectives 


This chapter describes the laboratory equipment and experi- 


mental procedures. The laboratory experiments were performed in both a 


straight and a meandering flume at the University of Alberta Civil 
Engineering Graduate Hydraulics Laboratory. The experiments had two 
primary objectives: 


1. To evaluate the effect of an ice-cover on the diffusion 


and mixing phenomena, and 


2. To provide information on the transverse diffusion and 
mixing characteristics in straight and meandering ice- 


covered channel flows. 


4.2 Experimental Equipment 


4.2.1 Straight Flume 


A general view of the straight flume is shown in the photo- 
graph of Figure 4.1. The flume was rectangular in cross-section, 
approximately 60 ft. long, 4.0 ft. wide, 2.0 ft. deep and had an ad- 
justable slope. The side walls were of plexiglass and the floor was 


plywood covered with a smooth layer of fiber-glass resin. 


Water was pumped from a storage sump into the flume head tank 


through a 12 inch overhead pipe fitted with valves and a recording 
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FIGURE 4.1 GENERAL VIEW OF STRAIGHT FLUME 
FLOATING ICE-COVER. 
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8 inch Foxboro magnetic flowmeter. Downstream of the flume the water 
flowed into a small channel and back into the sump. A 6 ft. long 
floating horizontal board was installed at the upstream end of the 
flume to smooth the flow. Depths of flow were measured manually along 


the centerline of the flume with point gauges. 


4.2.2 Meandering Flume 


The meandering flume was 62.5 ft. long, 2.5 ft. wide and 1.0 
ft. deep. Figure 4.2 is a schematic diagram of the flume. This flume 
had two identical but reversed 180° circular bends with a centerline 
radius of 9 ft. Because of insufficient laboratory space it was not 
possible to join the two bends with a straight section as naturally 
occurs in most rivers. The flume width of 2.5 ft. enabled width to 
depth ratios of up to 20 to be obtained. In this respect, the experi- 
ments in this flume were performed under conditions which approximated 


those of a natural stream. 


A 5 ft. long straight reach connected the first bend to the 
head tank. The flume bed was made with 3/4 inch plywood and the side- 
walls with 1/8 inch thick aluminium sheets. Joints were sealed with a 


silicon sealant and the entire flume painted. 


This flume rested on ten wooden supports, each fitted with a 
pair of 3/4 inch bolts for slope adjustment. It shared the same pump- 
ing, flowmetering and water storage facilities with the straight flume. 


At the inlet, water was discharged vertically downwards into the inlet 
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FIGURE 4.2a PLAN VIEW OF MEANDERING FLUME 


FIGURE 4.2b PHOTOGRAPH OF MEANDERING FLUME, 
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tank through an 8 inch pipe. A set of vertical screens was installed 


in the inlet tank to smooth the flow into the flume. 


In all experiments, flume water was recirculated during 
velocity measurements but run to waste during tracer concentration 
measurements. The sump water was replenished continuously from the 
city mains during the concentration measurements. Because discharge 
from the city mains was limited it was not possible to perform tests 


with large flow depths. 


4.2.3 Bed Roughness and Ice-Cover 


In general, the bottom resistance to flow in an ice-covered 
stream is different from that of the ice-cover. To simulate this 
condition in the laboratory, expanded sheet metal, 1/16 inch thick was 
used as bottom roughness in both the straight and meandering flumes 
(see Figure 4.3) for some experiments. The ice cover was simulated 
with painted 1/2 inch plywood boards stiffened along the edges with 


2" x 2" wooden strips to reduce warping (see Figures 4.1 and 4.4). 


At the connection between these plywood boards a horizontal 
spacing of approximately 1.0 inch was left to facilitate velocity and 
concentration measurements. The weight of the "ice-cover" was 2.0 lbs. 
per he It floated free on the water in both flumes. The "ice-cover” 
inevitably warped in spite of the wooden strips used to stiffen it. 
Experiments were therefore performed with flow depths greater than 


0.10 ft. to minimise the effects of the warping of the "ice-cover" on 


the flow. 


—_ 


ié 


LS 7 , 1 is q ott ‘s 
bolissent sev exovrod Inobazev io ten A -egte 


a? an 


se alt sank wold, add) e00me 


gaits bats (ine ciSSt.,2ay Keohane eens ie 
soo Soom aokaud siesw os oie wea 
aio wort ylavavmigaos ama nated Baw vga. . 
stxadoeth. sevanad . sadnans tienen norimninscnad Male k = 
atess miniteg o2 sitieeog Joo asw, o2 bestaht eae ome' 


. fi tee ie | 
barevoo-s4F nw At wal} og sonndbhesy hid oats see Le 
gids sinluwte of Tevoo-eak ots he an eh ages 

eew xobdavitone OR sit cutout: ‘5 
2en0l3 gabsntacay bs: a | 
beintumbs aie TINGD" ak BAe — cme mt 
Rebate abghe ant sire begets nbasad tie 
te a ie 


Letnasktion: ® gr imide naa 
baa eotsenay tra fer) oe. RS em, dant os 
-2dL 0.5 gape “sous ait oe te _axtow itt 
“yeveostiatt batt sll deli! ode , 
at agathae 03 saith ears cin lt 19 cate « we bagtev cidaateoat 2, 


lant, 


<i) x q 


ties aSohorg edsgeb wold Saad lee oro2ereds snew eaneatzeqes 1 J 
to to “ieves Heo”. ots te amie of 8 Ca ay — oF it ane 


By 


FIGURE 4.3 FLUME BED ROUGHNESS-EXPANDED SHEET METAL 
T/l6in. thick. 


FIGURE 4.4 PARTIAL VIEW OF MEANDERING FLUME SHOWING ICE-COVER 
(1/2 in. plywood boards) IN POSITION, 
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a a Velocity Measuring System 


4.2.4.1 Straight Flume: Velocities were measured with a 
1/8 inch diameter pitot-static tube. The difference between dynamic 
and static pressures was measured with a variable reluctance differen- 
tial pressure transducer with a 0.5 psid full-scale diaphragm (Pace 


model P1D). 


This pressure transducer was connected to a multi-channel 
carrier demodulator so that its output could be scanned using a Data 
Acquisition System. The transducer was calibrated by applying static 
differential heads of zero to 1 inch of water to the pitot-static tube 
and adjusting the demodulator for a linear output of zero to 10 volts. 
This calibration was checked before and after each experiment. The 
differential pressures obtained from the transducer were converted to 


velocities, without a turbulence correction, and then averaged. 


Because of the relatively slow sample interval of the paper 
punch system, each channel of the demodulator was provided with variable 
damping circuits to reduce the frequency response of the transducer 


output to a level suitable for input to the data acquisition system. 


4.2.4.2 Meandering Flume: Water velocities in the meander- 
ing flume were determined by a three-tube yaw probe capable of giving 
the magnitude and the direction of the velocity vector component in a 
horizontal plane. The principle of the probe is described later. The 


probe was similar to that used by Rajaratnam and Muralidhar {196732 0ne 
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was made of 3 mm 0.D./1.8 mm I.D. stainless steel tubing. At the nose 
of the probe the side tubes were chamfered at 45° as shown in Figure 


4.5. 


Differences in dynamic pressure between the two outer tubes 
and the central tube were measured with two (Pace model P90D) differen- 
tial pressure transducers equipped with a 0.03 psid full-scale dia- 
phragms. The calibration procedure for these transducers was identical 


to that for the straight flume experiments. 
4.2.5 The Yaw Probe 


Let q_ be the velocity vector at a point in a two-dimensional 
turbulent open-channel flow and w be the angle of deviation from a 
fixed reference direction. If the yaw probe is positioned in the refer- 
ence direction, the total heads that will be indicated by the three 


tubes can be expressed as: 


H = H +k, (4.1) 


where the number 2 stands for the centre tube, 3 the tube on the side 
of the velocity vector and 1 the remaining tube; Hy is the static 


head; Kp Ky and K3 are calibration coefficients and functions of 
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only wW for a given probe if viscous effects, velocities normal to 
the plane of the probe and other minor correction factors are regarded 


as negligible. 


The yaw probe was calibrated in the potential core of the 
plane wall jet with zero pressure gradient, produced by a deeply sub- 
merged sluice gate in a rectangular channel. The velocity in the 
potential core of the jet was determined by a pitot-static tube and 
the angle of yaw was. varied from -50 to +50 degrees during the cali- . 
bration. The calibration was performed for a velocity of 3.2 ft/s 


only. Care was taken to ensure that the plane of the tubes was hori- 


zontal. 
From Equation 4.1 a fourth coefficient Ky can be defined 
as: 
(Daa er ee 
De ee ae oo Kea (4.2) 
2 a 2 2. 


K, is again dependent only on y. The calibration curve given by 


4 
Rajaratnam and Muralidhar (1967), obtained for velocities ranging from 
1.47 to 10.39 ft/s, was found to describe the data very well and was, 


therefore, adopted. This curve is shown in Figure 4.6. The negative 


half of the calibration curve is given by symmetry. 


It is noted that the yaw probe was calibrated in a uniform 


velocity field away from boundary or wall effects, whereas the condi- 


tions during measurements were quite different. Near the flow bounda- 
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ries viscous effects coupled with significant velocity gradients, would 


inevitably introduce some errors. It is possible to minimize the effect 


of velocity gradients by using a smaller size probe. However, such a 
small probe would have had a very slow response. Rajaratnam and Mura- 
lidhar (1967) conducted some measurements with the yaw probe in a wall 
jet and found that the uniform flow calibration curve could be used 


for some shear flows without serious error. 


The procedure used for evaluating the velocity components 
u and w is outlined below. Ky, is first found from Equation 4.2 


and the corresponding value of  ~ is read from Figure 4.6. From 


Equation 4.1 it can be shown that: 


Ze (4) i=, ) 
pian Ob aie Saree | (3a) 
Poa aL 
or: 
= J 4.3b 
q Ks 2gAH (4.3b) 
Z 
where: kK. *: = and AH = H, - H 
5 4 me 2 A. 
> ie 
Ks is assumed to be a function of only. A plot of K, versus Ww 


is shown in Figure 4.7. Using the known value of ADs ng Ks is found 
from Figure 4.7 and the magnitude of the velocity vector, q , is com- 


puted from Equation 4.3b. The two velocity components are given by: 
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4.2.6 Tracer and Injection System 


The general procedure for determining the rates of transverse 
diffusion and mixing experimentally is: 


1. To introduce a tracer continuously at a steady rate into 


the flow at some point, and 


2. Measure the tracer concentration at various cross sections 


downstream of the source. 
The three kinds of tracers commonly used in water are salt, fluorescent 
dye and radioactive isotopes. Standard methods are available for mea- 


suring the concentration of all three. 


The tracer used was Rhodamine WT fluorescent dye, a 20% 


solution having a specific gravity of approximately 1.2. Concentrations 


of Rhodamine WT as low as 0.1 ppb can be detected with a fluorometer 


by very simple procedures. 


The tracer injection system was the same for both flumes. A 
sketch of the system is shown in Figure 4.8. The tracer was always 
diluted to less than 100 ppm with dechlorinated tap water and the solu- 
tion stored in a 30-gallon Mariotte tank. Dechlorination was necessary 
because Rhodamine WT at low concentrations is found to react with 
residual chlorine in tap water. The dechlorination was achieved by 


adding about 3 ppm Sodium Thiosulphate to the water. 
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FIGURE 4.8 SKETCH OF TRACER INJECTION SYSTEM 
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A precision rotameter capable of measuring very low flows 
indicated the tracer flow rate. The calibration curve for the rota- 
meter is shown in Figure 4.9. The tracer was released into the flow 
at a continuous and steady rate through a 1/8 inch stainless steel 
tube. In all experiments, the tracer was introduced at mid-depth paral- 


lel to the flow and at a velocity which approximated that of the water. 


4.2./ Concentration Measuring System 


4.2.7.1 Straight Flume: The sampling probe was a 1/8 inch 
diameter stainless steel tube fitted into a holder mounted on an alu- 
minium angle which was placed across the flume. The latter was provided 


with a scale for reading the lateral position of the sampling probe. 


Samples were obtained by siphoning the flume water into 60 cc 
bottles and then analysed for concentration with a fluorometer (Turner 
Model 111). The sampling velocity was adjusted to approximate that of 


the mean flow velocity. 


4.2.7.2 Meandering Flume: The first few tracer concentra- 
tion measurements in the meandering flume were conducted utilising a 
system identical to that employed for the straight flume experiments. 
This eseten was later replaced by a sampling rake having 24 sampling 
probes made of 1/8 inch soft brass tubes. The sampling rake is shown 
in Figure 4.10. This rake enabled simultaneous sampling at 24 points 


at 0.1 ft intervals across the flume. 
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FIGURE 4.10 SAMPLING RAKE USED FOR TRACER MEASUREMENTS. 


FIGURE 4,11 PHOTOGRAPH OF DATA ACQUISITION SYSTEM, 
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4.2.8 Coordinate Positioner 


The electronic coordinate positioner was the same as that 
used by Hollingshead (1972) and is basically a y-z "plotter" that 
positions a probe a any desired location within a cross-section. The 
DC voltages representing the coordinates were first zeroed with refer- 
ence to the flume side and bed as datum planes and then fed into the 


data acquisition system. 


An instrument carriage for the velocity probe was mounted on 
both the straight and meandering flumes, with the carriage being manu- 
ally positioned at any desired cross-section (x-coordinate) and the 
velocity probe electronically positioned within the cross-section. 


4.2.9 Data Acquisition System for Velocity 


Measurements 


The Data Acquisition System consisted of the following ele- 


1. 20 Channel Scanner (VIDAR 604). 

2. Integrating Digital Voltmeter (VIDAR 500). 

3. Data Coupler (Control Equipment Corporation Model 310). 

4. Scan Counter. 

5. Teleprinter (Teletype Model ASR 33). 
A photograph of this equipment is shown in Figure 4.11 and is repre- 
sented schematically in Figure 4.12. All inputs for the data acquisi- 


tion system were in the range zero to ten volts and each was assigned 
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to a particular channel of the Scanner. The Scanner sequentially 
sampled selected channels and routed the analogue signals to the Digital 
Voltmeter where they were digitized and converted to Binary Coded 
Decimal (BCD) form. The Data Coupler formated the BCD signal and direc-- 
ted it to the Teleprinter. The Teleprinter then punched the data on 
papertape in American Standard Code for Information Interchange (ASCII) 
and provided a decimal printout on paper. The desired number of scan 
cycles and the scan rate were selected on the Scan Counter and the 
Scanner respectively. The upper limit for the scan rate was one read- 
ing per second. This was determined by the permissible teleprinter 


speed of 10 characters per second. 


4.3 Experimental Procedure 


4.3.1 Identification Code for Flumes and Experiments 


The identification code for the flumes consists of two or 
three letter combinations. The first letter is either S or M, 
S representing the straight flume and M the meandering flume. The 
second letter in the code indicates the type of bed roughness. R 
denotes a flume bottom roughened with galvanised sheet metal and S&S 
a css bottom. The third letter, when present, is always I and de- 
notes the presence of an ice-cover. Thus the code MRI refers to an 


experiment conducted in the meandering flume with a roughened flume 


bed and an ice-cover. 
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All experiments were grouped in "series", each series being 
made up of a number of "Runs". The first digit of an experiment number 
refers to the series, the second and third digits to the run within the 
series. For example. run 304 stands for experiment 4 in the 300 series. 
Series 100 and 200 were conducted in the straight flume, series 300 
and 400 in the meandering flume. Classification of the flumes and 


experiments is summarised in Table 4.1. 


4.3.2 Description of a Typical Experiment 


4.3.2.1 Straight Flume: The first step in all experiments 
was to establish uniform flow conditions by setting the flume slope and 
adjusting the discharge and tailgate until the desired normal depth 
was achieved. Uniform flow conditions were assumed to prevail when 
depths measured along the flume centerline were within + 0.001 ft. of 


the average depth. 


Velocities were measured at seven verticals located at z/w 
= 0,0; 3.0.25.) 920.3759 and, £40.4375—g-and. LOgpoints in ‘the vertical 
at only one cross-section. This cross-section was located approximate- 


ly half-way along the flume length. 


The velocity distributions were determined as follows: 


1. The y and z coordinates (coordinate positioner) and the 


pressure transducer were zeroed and calibrated. 
2. The scan rate and the number of scan cycles were selected. 


3. The velocity probe was positioned at the desired point and 


the counter started. 
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4, Step 3 was repeated for all selected points within the 


cross-section. 


The pressure transducer calibration was checked after each run. 


The concentration measurements in the straight flume were 
made after velocity measurements were completed. In all the experi- 
ments, four measuring cross-sections were used. The concentration 
measurements were made at three levels; y/D = 0.25, 0.50 and 0.75. 
Sampling began at the cross-section nearest the tracer source and pro- 
ceeded downstream until all cross-sections were covered. The sampling 
period at each point was about 25 seconds. The total time required 
for sampling at the last section where the tracer cloud width was 


greatest was approximately 90 minutes. 


Tracer measurements were carried out for both centerline and 
side (3 inches from side-wall) injections under identical flow condi- 


tions (for example, runs 101 and 102). 


4.3.2.2 Meandering Flume: Again the tests were started 
with the establishment of uniform flow in the flume. For the open- 
channel tests uniform flow conditions were assumed to have been 
attained when the flow depths on the centerline were within + 0.002 ft. 
of the overall average depth. The cross-sectional average flow depth 
for ice-cover flows was obtained by taking the average of the depths 
at the left bank, right bank and channel centerline. Due to the 


warping of the plywood "ice-cover", deviations of the cross-sectional 
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average depth from the overall average flow depth could only be kept 


within + 0,006 £c. 


Velocities were determined at 5 verticals located at z/w = 
0, +0.20 and + 0.40, and 10 points in the vertical for the sections 
indicated in Figure 4.13. The tracer concentration measurements were 
made at three levels; y/D = 0.25, 0.50 and 0.75 for 8 to 11 sections. 
Three tracer source locations were used. These were at the channel 
eee tine and at 0.05 ft. measured laterally from both the left and 


right banks. 
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CHAPTER 5 


PRESENTATION OF VELOCITY MEASUREMENTS 


5.1 Straight Flume 


5.1.1 Hydraulic Data 


The hydraulic conditions for all test runs are summarised in 
Table 5-1. Column 3 lists the uniform flow depth, d, as measured 
along the channel centerline. The reference level for d, was taken 
at one-half the thickness of the galvanised sheet metal (bed roughness) 
above the flume bed. The average shear velocity in Column 8 was eva- 
luated from u, = VeRS , where R is the hydraulic radius; Se is 
the slope of the energy line and g is gravitational acceleration. 
The average velocity, Uy » was calculated by dividing the metered dis- 


charge, Q , by the cross-sectional area of flow. 


5.1.2 Velocity Measurements 


Velocity profile measurements were made at only one cross- 
section for each test run in order to check the metered discharge and 
to ascertain that the flow was two-dimensional. Velocities were 
determined at seven points laterally and ten points vertically. The 
discharge calculated from the velocity measurements was within -6% of 
the metered discharge. The integration of the velocity distribution 


near the flume walls contributed substantially to the discrepancy in 
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the metered and calculated discharges because of the assumption of 
linearity between any two adjacent velocity profiles. Also the flow- 
meter could be read only to the nearest 0.01 cfs which could represent 


an average error of about + 3% for the discharge range used. 


Secondary circulation was observed in all flows and within 
the central 50% of the channel variations in the depth-averaged velo- 
cities of up to 12% were recorded. For the open-channel experiments 
the velocity profiles were found to be logarithmic; Figure 5.1 gives 


the velocity distribution for run 106. 


The normalized mean velocity profiles near the center of the 
channel for all the test runs with ice cover are shown in Figure 5.2. 
The value of uF x in each run had to be taken from a plot of the 
velocity profile and this probably contributed to the small discrepan- 


Oo 
cies in the u/u ax versus y/D plots. 


The mean velocity distribution is distinctly asymmetric with 
the point of maximum velocity located nearer the ice-cover which is the 
relatively smoother boundary. Velocity contours for run 203 are plotted 
in Figure 5.3. They show that the flow was not entirely two dimensional 


but is more concentrated near both side-walls. 


From the logarithmic plots of the velocity profiles shown in 
Figure 5.4 it is evident that there is a region close to the rough wall 
where the velocity profile conforms with the logarithmic law. It can, 


therefore, be assumed that in this region the mean flow is unaffected 
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by the ice-cover. However, farther from the rough bed the velocity 
distribution is seen to deviate from a logarithmic relationship. The 
normalized distance from the bed, y/D » at which this deviation starts 
is about 0.4. Figure 5.5 gives the mean velocity distribution near 
the ice-cover. Although insufficient velocity measurements were made 
close to the ice-cover the data points seem to follow a logarithmic 


relationship within a narrow zone. This zone is relatively smaller 


than that near the rough bed and extends to about y/D = 0.2 , y being 


the distance measured downwards from the ice-cover. 


From the above it is evident that near both the ice-cover and 
channel bed the mean velocity profile is unaffected by the opposite 
boundary. The effect of the rough wall extends farther into the flow 
than that of the smooth wall. Between these two regions is a transi- 
tion zone wherein the velocity profile depends on the nature of both 
boundaries. These observations are in accord with those of Hanjalic 
and Launder (1968) who carried out their measurements in a wide rectan- 


gular conduit with one of the principal walls roughened. 


Because the number of point velocity measurements in the 
vertical were insufficient it was not possible to accurately determine 
the location of the point of maximum velocity and to evaluate whether 
the Reynolds Number, Ry , has any effect on its position. However, 
Hanjalic and Launder observed that increasing the Reynolds number moves 


the maximum velocity closer to the smoother boundary. They also found 


a similar shift in the position of the zero-shear plane which was always 
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nearer the smoother boundary than the point of maximum velocity. 


5.2 “Characteristics of Flow in Meandering Channels 


5. 2. 1. wGeneral 


Flow in meandering channels differs from that in straight 
channels because the presence of centrifugal forces results in the for- 
mation of transverse pressure gradients and the associated lateral 
pial aeton in the plane of a cross-section. In general, the flow 
structure changes in the downstream direction due to the growth, decay 


and reversal of this circulation or spiral motion. 


In order to simulate flow conditions in natural streams the 
present experiments were conducted for W/d, ratios greater than 10.0. 
The Froude Number was also kept low, the highest being 0.53 and the 
Reynolds Number was greater than 10° in all tests so that they were 
within the turbulent regime. A summary of the heirene data is pro- 
vided in Table 5-2. For each test run, the discharge through any cross- 
section used for velocity and tracer measurements was calculated from 
the velocity measurements. This calculated discharge was found to be 
within + 8% of the metered discharge. The discrepancy is attributed 
to the small number of vertical velocity profiles taken, measurement 
errors and the assumption of linearity between adjoining velocity pro- 
files. Also, as stated previously, the magnetic flowmeter could only 


be read to the nearest 0.01 cfs which could cause an average error of 
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about + 2.5% in the discharge range used here. 


5.2.2 Velocity Measurements in Open Meandering Flume 


Mixing in open and ice-covered channels depends greatly on 
the velocity distribution. The spiral motions generated by bends do 
not only modify the velocity distribution but also augment lateral 
mixing. Detailed measurements of the velocity field are, therefore, 


necessary for an understanding of the mixing process. 


5.2.2.1 Longitudinal Velocity Component: Typical plots of 


longitudinal velocity measurements using the data from run 301 is pre- 
sented in Figures 5.6a, b and c. The corresponding measurements for 
runs 303 and 307 are shown in Figures A.1 and A.2 of Appendix A. A 
study of the measurements reveals that close to both the inner and 
outer banks, the maximum longitudinal velocity does not in general 
occur at the surface but is displaced downwards. This is caused by the 
spiral motion which transfers low momentum fluid from the bottom 
layers to the water surface near the inside bank, while near the outer 
bank it supplies high-momentum fluid from the surface downwards. Also, 
the formation of a circulation cell with an opposite sense of rotation 
to the main spiral motion at the upper corner of either the outside or 
inside of the bend could shift the maximum longitudinal velocity down- 
wards. The reduction in velocity near the free surface at the outer 


bank appears to be smaller than that near the inner bank. 


The transverse distribution of the depth-averaged longitudinal 
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velocity, u , are aie As Figures 5./a, pyre. It is evident from 
these figures that the higher velocities occur near the inside of the 
first bend at the beginning of the bend. This observation is in agree- 
ment with potential flow theory for flow around a bend which predicts 

a free vortex type of motion with u ~ l/r. The high velocity zone 
starts to shift towards the channel centerline after angle, 01 » from 


the beginning of the bend exceeds about 1/6. 


Downstream from oT = 1/6 the zone of maximum velocity gra- 
dually moves away Br Ae ie inside bank and at oy = 27/3 Sit is near 
the outer bank. On emerging from the first bend the zone of high 
velocity flow is located near the inside of the second bend. This pro- 
duces a marked non-uniformity in the lateral distribution of u which 
persists in the second bend until oP) = 1/3 , where oP) is the angle 
in the second bend. From this section downstream the high velocity 
filament again moves across to the outside bend. If this series of 
bends were continued, odd-numbered bends for example, will exhibit 
similar depth-averaged longitudinal velocity profiles at corresponding 


sections. 


The experiments also indicate that the displacement of the 
high pelociey flow region from the inside to the outside of a bend is 
accomplished over a relatively shorter distance for flows over rough 
beds compared with that over the smooth bed. The cause of this is not 


very clear but is believed to be related to the non-uniform distribu- 


tion of bed resistance. 
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A parameter which can be used to describe the non-uniform dis- 
tribution of velocity is the dimensionless mean square velocity devia- 
tion, 6 , defined as: 

+W/2 
(u - U,)*az 
1 -w/2 
Hi 2 
where W is the channel width and Uy is mean velocity. The longi- 
tudinal variation of 6 for all tests is shown in Figure 5.8. The 
contributions to § are due only to velocity variations in the lateral 
direction. It is seen from Figure 5.8 that high values of 6 occur 


between 0, = 27/3 and 6, = 1/3 and in run 307, where the channel bed 


3 
was smooth, 4 retains a high value even at 0. = 1/2 because the 
high velocity zone was still located within the inside half of the 


bend at that section. The degree of non-uniformity in the velocity 


field was observed to be relatively low within the first bend. 


5.2.2.2 Transverse Velocity Components: Measurements of 


the transverse velocity component from the data of run 301 are shown 
plotted in Figures 5.9a, b and c. The corresponding measurements for 
runs 304 and 307 are given in Figures A.3 and A.4 of Appendix A. The 
plots clearly show the growth and decay of transverse circulation or 
spiral motion in the two bends. Following Yen (1972), the strength of 


the spiral motion, S , can be defined as: 
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S therefore denotes the sum of the inward and outward unit lateral dis- 
charges and is very convenient for computational purposes. The longi- 
tudinal variation of S is presented in Figure 5.10, where the solid 
lines describe the centerline variation of S and the plotted points 


represent the variation of S within the central 40% of the flume. 


Since the sense of rotation of the spiral motion reverses in 
the second bend the spiral motion strength, S , in this bend should be 
= ee differently from that in the first bend. However, because 
the number of lateral point velocity measurements was small it was not 
possible to account properly for the residual spiral motion from the 
first bend in the second bend. Therefore, in Figure 5.10 S is plotted 
with the same sign for both bends. The stream function cannot be used 


as an alternative parameter because the flows are three-dimensional. 


The strength of the spiral motion varied laterally in all 
cross-sections as illustrated in Figure 5.11 for run 407. When this 
figure is compared with Figure 5.14c it is seen that the shift of the 
region of highest spiral motion strength from one bank to the other 
generally follows a pattern which is similar to that of the maximum 


velocity filament. 


It is evident, from Figure 5.7 and from the fact that there 
was neither a section with zero lateral discharge nor any flow para- 
meters independent of x , that the fully developed bend flow condition 


was not attained in any of the test runs. 
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FIGURE 5.10 LONGITUDINAL VARIATION OF THE STRENGTH OF 
SPIRAL MOTION NEAR CHANNEL CENTERLINE FOR 
OPEN CHANNEL EXPERIMENTS 
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All the test runs indicate the presence of transverse circu- 
lation at the beginning of the first bend and the spiral motion develops 
relatively faster in the rough bed experiments, runs 301 and 304, 
compared with the smooth bed experiment, run 307. It attained its peak 
near 85 = 1/3 for the rough bed experiments and at 05 = 27/3 for 
the smooth bed test. Further downstream, but still in the first bend, 
the ray motion gradually starts to weaken. Yet at 9, =T or 0, 
= 0 it still retained considerable intensity. Minimum values of §S , 


the strength of the spiral motion, were recorded between 0. = 0 and 


oP = 1/6. 


Downstream of 0, = 0 the helical motion due to the second 
bend begins to grow and occupy an increasing area of the flow cross- 
section (see Figure 5.9 and Figures A.3 and A.4 of Appendix A). At the 
same time it displaces the residual spiral motion from the first bend 
towards the outside of the bend removing it completely between 0, = 1/3 


and 6, = 27/3. 


2 


Thetteoul ta tend to show that the decay of residual circula- 
tion from the first bend is controlled by the roughness and width-to- 
depth ratio. The residual spiral motion persisted over a longer reach 
for en smooth bed flow than for the rough bed flows. It was also 
observed that the residual circulation decays relatively faster for 
large W/d, ratios. These observations are in agreement with the find- 
ings of Yen (1965). The helical motion due to the second bend did not 


reach its peak strength until after the residual spiral motion had 


99 


~ \ : 
: nai i 
~yotts aetevensyd 4 vomeeeq ode atertbal acer oy 
ecoloveb tolsos faxtqe nda Bis bed gazl? ot ~e 
aie: 
oO baw. [Of enw eanemiatecil bead y007 arid ra Kis he 
; : % 
$294 edt bonkatite ‘7 “Oe sine. sremiteqxe bad soon wt i 


oe t\ni « 6 de hoe esesetkragne bad dgyor, at = tes . yw 


Baed test? nild a8 cies sud ceusrtinwob zonjxot ead | ve 
8 a0) n= Boe sat’ -raskn git “9 sityese iasbert, ‘wo som | tea 


Jo Renta ‘wemtnth ee keene k biderebisnoo bantetos Litse | 
hon «20 = <P ngawiod bwlsaqied' ane notaee ‘tiostge ot to 7 


os9e of? of sub cotrom Deotiod ane 0 . a aoe 
~2¢019 ole of Yo tad as has Se ied gqe000 bas vores 2 


1538 . Ch-wkbnsgqaiit®, Babee £14 Soaegt bre ee omgrt ~~ 
. booed 2078, oft most nabsea Peabye foubtems 383 connie ats 
e\s m8) meswaad eee as aes bisd add to nee ot 


at mi ee 


v 


— 


“siyouls Lautner 20 yodeb aid asia wore 03 baat ised: PPS 
~os-~4 bie bee naclligess only 4a. palttexanes at baad seat? ans" 7 
dsaexy zagnol a teva beagteved aobsaiefartge | feobkaex od ea | 

coals eew 37 awa’ hed dguot. eis 302 oads ro bod aa 

zt 299en) Fisvidesor eypoeh nobielworts isutitess 62 det bovreede 
. bab? att dziw tooreatgeh nt ESTs enoidaviseda ovedt seeksex jw ogral 
~ Hon bkb bred bnosaa +l ay sub: polzam Geotied sit eee) nox to oe 


oar, Gh tiptoe. Scag Ibe oat vests ihaaw coe deen eat tone 


fs a 
t - 
J i 
hi 

5 


eee 
at 7 iy. i 
va ‘marl te 


mete 7 eh 
a a a 
D> ae) 
7 


100 


decayed completely. It then decays towards the bend exit because of 


the presence of the tail-gate. 


Maximum observed transverse velocities near the bed were 
slightly larger in the first bend than in the second bend. On the other | 
hand the transverse velocities near the free surface were approximately 
the same for both bends except in run 307 where higher transverse sur- 
face velocities were recorded in the second bend. Most vertical profiles 
indicated the reduction in the transverse velocity near the bed begin- 
ning further from the bed for the rough bed flows (runs 301 and 304) 
than for the smooth bed flow (run 307). As could be anticipated, the 
transverse velocities are smaller near the side-walls than within the 
central portions of the flow. Near the walls the vertical velocity com- 
ponents become important and continuity considerations require that the 


w-component decrease. 


5.2.2.3 Angle of Deviation of the Velocity Vector: The 


variation of the horizontal angle of deviation of the velocity vector, 

Ww , in the vertical and lateral directions is presented in Figure 5.12 
and Figures A.5 and A.6 of Appendix A. The angle w can be regarded as 
an indicator or measure of the spiral motion since yp = Func (anus The 
en inward angle occurred near the channel bed and the maximum out- 
ward angle near the free surface. In general Ys! > Yea for most 
verticals. The zero angle of deviation in a vertical was located in the 


lower half of the depth for practically all verticals. 
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According to Yen (1972) and Rozovskii (1963) tan (dy me 
d & 
* 
varies approximately with the product i ao. (where K is 
Crakk 


von Karman's constant). Rozovskii recommends that for both smooth and 
rough channel beds tan ~ near the bed be approximated by (10-12) d,/r 
for fully developed bed flow. The observed maximum values of ee 


Oo 


are compared with Rozovskii's approximation in Table 5-3. 


TABLE 5-3 COMPARISON OF OBSERVED AND COMPUTED Tee, VALUES 


d 


Run i oe as. Experimental 
Number OAs OR eoten | [| (10-12) z_! Remarks 


301 1/39 24a em oy Le Rough bed 


304 1/64.7 Suc. pO. 5 - Rough bed 


307 1/38.3 14.6° = 17.4° Smooth bed 


The values of Wes observed in this study are of the same order of 
magnitude as those given by Rozovskii's relation. However, the observed 
~ values contain measurement errors due to inaccuracy of the yaw probe 
in steep velocity gradients. Also it is noted that the width to depth 
ratio does not appear in Rozovskii's approximation for Me 3 and yet 
experimental evidence (Rozovskii, 1963) indicates that it has a signifi- 
cant influence on the development and decay of the spiral motion and 
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5.2.3 Velocity Measurements in Ice-Covered Meandering Flume 


One characteristic difference between flows in ice-covered 
and Shel EE bends is the formation of a double transverse circula- 
tion or spiral motion in ice-covered bends due to the presence now of 
velocity gradients at both the bed and the ice-cover. In this section 
some properties of flow under ice-cover conditions are described and 


compared with those observed for open water conditions. 


5.2.3.1 Longitudinal Velocity Component: The longitudinal 


velocity measurements are presented graphically in Figure 5.13a, b, c 
for run 404. The velocity data for other test runs in series 400 are 
given in Figures A.7 and A.8 of Appendix A. Near the outside of the 
bend there is transport of high momentum fluid from the central portions 
of the flow to the top and bottom layers. This is expected to produce 

a fairly uniform distribution of the longitudinal velocity over a sub- 
stantial portion of the depth on this side of the bend. The same effect 
on the longitudinal velocity distribution near the inside of the bend 

is anticipated by the transfer of low momentum fluid from the top and 
bottom layers to the central portion of the flow. The above observa- 
tions were found to be valid for many cross-sections especially near the 
freidachend for runs 401 and 407 but not for run 404. In this run, the 
transverse velocities were much less than the longitudinal velocities 
and consequently the transport of low momentum fluid by the spiral 
motion to the central portion of the flow did not cause any noticeable 


change in the longitudinal velocity profile near the inside of the bend. 
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FIGURE 5.13b VERTICAL DISTRIBUTIONS OF LONGITUDINAL VELOCITY 
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For the rough bed experiments, that is runs 401 and 404, the 
maximum velocity was generally located between 0.50 to 0.70 of the local 
depth from the channel bottom. With a smooth bed it was close to mid- 


depth (Figure A.8). 


The lateral distribution of depth-averaged longitudinal 
velocity along the two bends is given in Figure 5.l4a, b, c. The Rien 
velocity portion of the flow is again initially located near. the fnatde 
8, = 1/6. On emerging from 


the first bend the fluid near the inside of the second bend is further 


. 


bend but starts to depart from it after 


accelerated and this results in a pronounced non-uniformity in the 
lateral distribution of u Beha at to that observed for the open chan- 
nel tests. Downstream of 0, = 1/3 the high velocity filament again 
gradually crosses over to the outer bend. The displacement of the high 
velocity area from one side of the bend to the other was observed to be 
relatively more gradual with an ice-cover than the corresponding open 


channel tests and is probably due to increased resistance to the flow. 


The variation of the dimensionless mean square velocity devia- 
tion, 6 , is presented in Figure 5.8. In general, the longitudinal 
variation of 6 was similar to that observed for the open channel tests 
Ee within the second bend. Also the average value of § for 
run 401 (rough bed) was higher than that for run 407 (smooth bed). Both 
runs had approximately the same W/d, and US: The same was observed 
for runs 301 and 307 in the open channel tests. This phenomenon is pro- 


bably caused by the combined effects of the spiral motion and the non- 
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uniform distribution of ice-cover and bed shear. The non-uniformity of 
the transverse depth-averaged longitudinal velocity profile increased 
with width-depth ratio for both ice-cover and open water experiments 


but the variation was smaller for ice-cover flows. 


5.2.3.2 Lateral Velocity Component: The measured local 


lateral velocity components are shown in Figure 5.d5a, b, ¢ for run 404 
and Figures A.9 and A.10 of Appendix A for runs 401 and 407. They 
illustrate clearly the growth and decay of the spiral motion. Two spi- 
ral motions with opposite sense of rotation were observed for most cross- 
sections. The bottom spiral motion generally occupied more than one- 
half the flow depth for the rough bed experiments. In the smooth bed 
experiments, run 407, each spiral motion was confined to approximately 
one-half the flow depth. The longitudinal variation of the strength of 
the spiral motion, S , is presented in Figure 5.16. The S-values for 
the ice-cover experiments include the effects due to the warping of the 
plywood used as the ice-cover. Substantial warping was noticed in some 
sections of the ice-cover. At 8, = 1/2 of run 404, for example, (see 
Figure 5.15) the greatest depth of flow was located on the inside in- 


stead of the outside bend as would have been expected. 


Transverse circulation was present at 0) = 0 and is obviously 
a manifestation of the lateral pressure gradient associated with the 
curvilinear nature of the flow. The spiral motions attained peak strength 
between oT = 1/6 and 85 = 1/3 in runs 401 and 404. However, for run 


407 (smooth bed) the spiral motions did not reach maximum strength within 
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0 0,1 0 et 0 Ort 0 Gl 0 O.l 0 OL Osz 
w in ft/sec. 
FIGURE 5.15a VERTICAL DISTRIBUTIONS OF TRANSVERSE VELOCITY, Runs 
404/5/6. SECTIONS 8,: O fT TT. 
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0.10 0 0.10 0 0.10 0 0.10 0 0.10 0 0.10 0.20 
w in ft/sec. 
FIGURE 5.15b VERTICAL DISTRIBUTIONS OF TRANSVERSE VELOCITY, Runs 
404/5/6. SECTIONS 9%, = 27/3 To $= Te. 
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0.10) 0,20 


w in ft/sec. 
FIGURE 5.15c VERTICAL DISTRIBUTIONS OF TRANSVERSE VELOCITY, \ Runs 
404/5/6. SECTIONS ®, = Ws TO a2w/3. 
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the first”bend until 6, = 27/3. The variation of S for run 404 was 


Z: 
rather erratic and is believed to have been caused by the warping of 


the ice-cover which had a marked effect on this run because of the small 


depth of flow. 


The spiral motions diminished downstream of the section of 
maximum S to a minimum value at Oo, eed Or 0 = 0 for all test runs. 
The gradual reduction in S within the middle of the first bend in run 
401 (Figure 5.16) is probably caused by increased resistance to flow 
because of the presence of the ice-cover while the further drop in S 
at the exit of the first bend, for all test runs, is due to the presence 
of the second bend. In contrast to the open channel tests the spiral 
motions due to the second bend were observed to have penetrated the 
entire flow width at 0, = 7/6. This actually is not surprising be- 
cause the average value of S within the first bend for the ice-cover 
experiments were about one-half the value for the corresponding open 
channel tests. Note that any two corresponding tests have the same 
experimental number in their respective series, e.g., runs 301 and 401. 
The new spiral motions in the second bend develop rapidly and peak be- 
tween oP = 1/6 and 8. an os They, however, decay near the bend exit 
because of the presence of a tail-gate at the end of this bend. By 
taking the average value of S for the two bends it is concluded that, 
for the W/d, ratios used in the tests, the strength of the spiral 


motion is greater if the width to depth ratio is smaller. This conclu- 


sion was also found to be valid for the open channel experiments. 
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Run 407 
0.02 ; W/d,= 10.2 


0.01 


0.01 


Run 401 
W/ d, = 17.0 


0.02 


Strength of the Spiral Motion, S ( ft2/sec) 


O T/6 1/3 Tr/2 57/6 T/6C/3 7/2 51/6 
2771/3 7,0 2717/3 


Section Number 


FIGURE 5.16 LONGITUDINAL VARIATION OF THE STRENGTH OF SPIRAL 
MOTION, S , NEAR CHANNEL CENTERLINE FOR ICE-COVER 
EXPERIMENTS 
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5.2.3.3 Horizontal Angle of Deviation of Velocity Vector: 
Measurements of the local horizontal angle of deviation of the velocity 
vector, ) , are given in Figures5.17 for run 404 and Figures A.11 and 
A.12 of Appendix A for runs 401 and 407. The largest angle of devia- 
tion occurred near the channel bed. For runs 401 and 407 the observed 
Ve oo | was 12° while that in run 404 was about 10°. These values are 
less than those observed for open channel tests and indicate that the 


strength of the spiral motion is comparatively less for ice-covered 


flows. 
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FIGURE 5.17a VERTICAL DISTRIBUTIONS OF HORIZONTAL ANGLE OF DEVIATION 
OF THE VELOCITY VECTOR, Y , Runs 404/5/6. SECTIONS 
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CHAPTER 6 


PRESENTATION OF CONCENTRATION MEASUREMENTS 
AND DISCUSSION OF RESULTS 


In this chapter the experimental results relating to tracer 
concentration distribution in the straight and meandering flumes are 


presented and discussed. 


6.1 Straight Flume 


6.1.1 Transverse Distribution of Tracer Concentration 


Typical transverse tracer concentration profiles at the four 
measuring sections are shown in Figures 6.la to d. The measured tra- 
cer concentration was normalized with the concentration for fully mixed 
conditions and no tracer losses, oe » and then plotted against the 
normalised transverse distance z/W for various x-values. With source 
at the a Sl centerline, the data indicate displacement of the cen- 
troid of the concentration profile from the channel centerline at some 
sections. No explanation can be offered for this because velocity 
measurements were made at only one section but it is speculated that 
this feature is probably caused by a tendency for the flow to meander 
along the channel. The shift in the centroid of the concentration 
distribution can also be seen for ice-covered flows and is more pro- 
nounced in these flows probably because of the warping of the plywood 


used to simulate the ice cover. 
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FIGURE 6.1 LATERAL DISTRIBUTIONS OF DEPTH-AVERAGED CONCENTRATION 
AT VARIOUS DISTANCES DOWNSTREAM OF SOURCE FOR 
CENTERLINE INJECTION 
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FIGURE 6.1 LATERAL DISTRIBUTIONS OF DEPTH-AVERAGED CONCEN- 
TRATION AT VARIOUS DISTANCES DOWNSTREAM OF 
SOURCE FOR SIDE INJECTION 
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When tracer is released near the side the peak depth-averaged 


concentration is observed to shift to the wall with longitudinal distance 


(see Figures 6.lc and d). This feature can be explained by the reflec- 
tion or image source technique (Sayre and Chang, 1968) in which the 
portion of the concentration distribution falling beyond the side is 
folded back and added to the portion falling inside the channel. It has 
been established that for uniform turbulent flow in a wide rectangular 
channel away from the side-wall boundary layers, the lateral concentra- 
tion distribution within a tracer cloud follows the Gaussian distribu- 
tion. This is confirmed in Figure 6.2 where the normalized tracer con- 


centrations are compared with Equation 3.18, that is: 


U,A E47 

« ols ; 
i = f (a) =: e 4k (3.18) 
C 

1m 

This equation may be rewritten as: 
b2 

C a Lee uk (2 Zz 
Ce ee eae ct he. (6.1) 
Cc = b 

lm ¢ 

m 


where is maximum depth-averaged concentration and b represents 
half the transverse distance between the verticals on which c = 0.5c 
In plotting these curves, transverse distances were measured from the 
vertical of maximum depth-averaged concentration. In Figure 6.2 the 
data points are plotted together with the Gaussian distributions ob- 


tained from Equation 6.1 using the measured values of iy and Er. It 
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FIGURE 6.2 GRAPH OF NORMALISED CONCENTRATION AGAINST NORMALISED 
DISTANCE: RUNS 106 AND 206 
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is evident that the fit is good for both open and ice-covered channels. 


6.1.2 Moments of the Transverse Tracer Flux Distribution 


6.1.2.1 Calculation of Moments: Lateral variation in the 
depth-averaged longitudinal velocity, a of up to 12% across the 
width of the tracer cloud was observed in all test runs. To allow for 
this variation of u the measured concentration values were discharge 
weighted so that the moments were actually calculated from the tracer 
flux (D uc) distribution instead of the usual concentration distribu- 
tion. The second moments about the point of tracer release, M, (x) A 


were calculated numerically by: 


N Vi N 
M, (x) = ty Ne ae MY (6.2) 
eZ (2 Jel tal oi 


in which Mi and Ms are the zeroth and second moments, respectively, 


of the tracer flux between any two successive data points about the 


point of tracer release. Mai and My are defined as follows: 


mee B Lain ECan) +THpecnp| Ge ti) (ors) 
and: Mos = D Venza CMgay) + Fapecang| 
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where nN and ns are the transverse distances of any two succes- 


igri 
sive data points from Zz the point of tracer release. The moments 


calculated by the above procedure are given in Table 6-1. 


6.1.2.2 Growth of the Second Moments with Longitudinal Dis- 


tance: Figures 6.3 and 6.4 show the growth of the second moments of 
the tracer flux with longitudinal distance downstream of the tracer 
source for centerline introduction of tracer. The plots indicate a 
linear growth rate of the second moments except close to the source. 
Since the centroid of the tracer flux distribution at various sections 
was found to be located very close to the channel centerline, the 


: . ; 2 
variance of the tracer flux distribution op, (x) was assumed equal to 


M, (x). 


The linear growth rate of of (x) extends nearly back to the 
source. This suggests that the lateral diffusion process converges to 
fre one-dimensional diffusion model fairly rapidly. The deviation from 
linearity close to the source is primarily due to the non-uniform 
vertical distribution of tracer concentration as shown, for example, in 
Figure 6.5 and to the fact that the tracer cloud might be too small for 


gradient-type diffusion to be truly applicable. 


From Equation 6.1 which was derived with the assumption of 
uniform velocity of flow across the channel, we can express the vari- 


; : Z 
ance of the tracer concentration distribution, ox) » as follows: 
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TABLE 6-1 NUMERICAL VALUES OF X, M, (X) AND A_ IN EQUATION 6.5 


CL refers to centerline injection. 
B refers to near-bank injection. 
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_ FIGURE 6.3 GROWTH OF VARIANCE OF TRACER FLUX WITH DISTANCE DOWN- 
STREAM OF SOURCE FOR OPEN CHANNEL TESTS 
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FIGURE 6.4 GROWTH OF VARIANCE OF TRACER FLUX WITH DISTANCE DOWN- 
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o (x) i= z + 6.4) 


We can conclude from Equation 6.4 and Equation 3.16 that b* o a, (x) 
and as b is a characteristic dimension of the tracer cloud width, 


that the cloud width grows parabolically with longitudinal distance. 


For tracer releases near the side-wall, the plots of M, (x) 
versus x given in Figure 6.6, do not clearly show a linear growth 
rate of M, (x) with x over the test reach. There is an initial 
reach over which M, (x) shows a fairly rapid growth rate, and this is 
probably due to the fact that the tracer was injected about 3 inches 
from the side and, therefore, it spreads laterally in two directions 
until no concentration gradients exist in one of these directions. As 
the tracer cloud spreads and occupies more of the uniform flow region, 
the growth rate of the second moments should approach the same rate as 
observed for centerline injection of tracer. It is not very clear 
from the data how quickly the linear growth rate of M, (x) is estab- 


lished in this case. 


6.1.3 The Transverse Exchange Coefficient 


6.1.3.1 Determined by the Method of Moments: It was shown 


in Section 3.3.1 that for tracer released from an arbitrary point with- 
in a cross-section, the average transverse exchange coefficient within 


a reach can be calculated using Equation 3.31, that is: 
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This equation can be rewritten as: 


dM, (x) 
= ee Nematic 
Ete Bio ade (6.5) 
W W 


2 2 
where M, (x) = | | uc ay | ucd 4] 


oa) 1 
Wo 2 

and A 4, cae [ fe en] 
Wy i} 


M, (x), A and x are ‘tabulated in Table 6.1. 


The results pertaining to centerline release of tracer are 
first considered. As expected A was essentially constant for a par- 
ticular test run with a variation not exceeding + 0.7% from the mean 
for the four measuring sections. Since E. and A are constant 
M, (xx) must grow linearly with x as observed in Section 6.1.2.2. 
Straight lines were, therefore, fitted to the M, (x) versus x plots 
and their slopes determined. Equation 6.5 was then used to calculate 


Ez. 


For near-bank release of tracer, A is no longer a constant 


but a function of x. The variation of A for runs 102, 104, 202 and 
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204 are shown in Figure 6.6 together with that of the second moments, 
M, (x). Equation 6.5 was used to determine average values of E. be- 
tween any two adjacent cross-sections. This average E, is plotted 

at the downstream section as indicated in Figure 6.7. The broken lines 
in Figure 6.7 are the E. values obtained for centerline introduction 
of tracer. The average E. between the tracer source and the first 
section of tracer measurement is lower than the average for each of 

the runs. Very close to the source only the small size eddies contri- 
bute to the difensaun of the tracer. As x increases and the tracer 
cloud expands, the average eddy size capable of making a contribution 
to the spreading of the tracer also increases resulting in a higher 
rate of diffusion. Since the largest eddy size is controlled by the 
flow dimensions and is partially dependent on the proximity of the side 
walls the rate of diffusion cannot increase indefinitely but should 
approach a limiting value. The smaller values of E. observed 
between the source and the first section of tracer measurement are, 


therefore, reasonable. 


The functional forms of A _ and M, (x) are not known and 
errors are inevitably introduced by assuming a linear variation of both 
variables between any two successive cross-sections to evaluate Eo: 
However, the technique used by Holley et al (1972) could be applied 


here by integrating Equation 6.5 once with respect to x to give: 


M, (x) = 2E., I(x) + constant (6.6) 
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where I(x) = x, It is to be noted that a plot of M, (x) versus 
I(x) must yield a straight line with a slope of 2E Some of the data 
for side release of tracer were analysed using Equation 6.6 but curves - 
were obtained when M, (x) was plotted against I(x). This could be 

caused by variable transverse diffusion coefficient near the side walls 


and the large distance between sections of tracer measurement. 


6.1.3 2; Determined by the Similarity Analysis: On the as- 


asumption of similarity, which is confirmed in Figure 6.2, the distri- 
bution of C,/C, or die~ was shown to be Gaussian in Section 3.1 
and the exchange or diffusion coefficient was derived as: 
bu 
me fe) 
Zz 4x1n2 


(3.20) 


It is to be noted that this equation does not apply in situations with 
side-wall effects. Equation 3.20 was also used to evaluate the average 
exchange coefficient over a reach of length x , where Us is now the 
mean velocity across the tracer cloud width. Only one measurement of 
the transverse tracer concentration profile is required to determine b 
in Equation 3.20. In this study the concentration profile obtained at 
the last measuring section was used to determine b and hence Eo 

The E. values so obtained together with those determined by the method 
of moments are listed in Table 6-2. The table also gives the percentage 
difference between the two E. values based on the E. values calcu- 


lated by the method of moments. The difference between the E. values 
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TABLE 6-2 COMPARISON OF E. VALUES AS DETERMINED BY 
THE METHOD OF MOMENTS AND SIMILARITY ANALYSIS 


Transverse Exchange 
Coefficient x 


E. ;: Bee ie Difference 


0.00092 . 00094 
0.00167 00174 
0.00134 .00137 
0.00142 .00149 
0.00050 00051 
-00075 . 00076 
-00074 0.00083 
.00067 0.00071 


as determined by Method of Moments 
as determined by Similarity 
Analysis 

E, (Similarity) ~ E, (Moments) 


E, (Moments) x 100% 
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was within 12%. In general, the differences were much smaller for 
open-water conditions than for ice-cover conditions. These differences 
can partly be traced to the assumption of constant uniform velocity | 
across the tracer cloud width. This simplification was not fully sa- 
tisfied in all the runs. Another cause of the differences is experi- 
mental errors. The average absolute difference between the E. values 
was about 4.5%. This indicates that Equation 3.20 is a satisfactory 


method for evaluating Eo 


6.1.4 Effect of a Smooth Ice-Cover on the Exchange 


Coefficient 


In order to evaluate the influence of an ice-cover on the 
diffusion process, experiments were conducted in both open and ice- 
covered channels. The results from the open-water experiments served 
as the basis for this evaluation. The flow depth and discharge were 
kept approximately the same for any two corresponding experiments (see 
Table 5-3). To facilitate comparison, any two corresponding experi- 
ments are given the same experiment number in their respective series. 


For example, run 201 corresponds to run 101. 


It is evident from the summary of data in Table 6-3 that the 
diffusion capacity of ice-covered flows, as measured by E. >» is sub- 
stantially reduced. An average reduction of approximately 50% was 
observed in all the experiments. The presence of an ice-cover reduces 


the significant size scale of the eddies and as a consequence there is 
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a reduction in diffusion potential. An obvious consequence of this is 
that the attenuation of peak concentration, which is indicative of the 
rate of spreading of the tracer cloud, is reduced in ice-covered flows 
compared with open-water flows. This also means that the channel 
reach required to achieve the same degree of mixing will be much longer 
for ice-covered flows than for the open-water case. These aspects of 


diffusion will be discussed later. 


6.1.5 Factors Affecting the Normalised Exchange 


Coefficient 


Turbulence in channel flows is generated by shear originating 
at the channel boundaries and this turbulence provides an efficient 
mechanism for mixing. Two features of turbulence relevant to diffusion 
are the turbulence intensity and the average size of the turbulence 
eddies. Laufer (1951) has oe that the r.m.s. value of the trans- 
verse velocity fluctuations is proportional to the shear velocity, u, , 
in a 2-dimensional flow in a wind tunnel. The size of eddies in a tur- 
bulent channel flow is limited by the average flow depth, d, , and the 
flow width, W. The above observations suggest normalising the exchange 
or diffusion coefficient with the product u,d,. The length scale, d, , 
prestdane regarded only as reflecting the average size of the dominant 


eddies producing the maximum mixing. 


Transverse velocities resulting from changes in channel geo- 


metry with longitudinal distance can affect E. » however, only rectan- 
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gular channels are considered here. 


Bulk properties of a uniform turbulent flow in a straight 
rectangular channel (e.g., the average transverse exchange or diffu- 
sion coefficient, E.. ,» friction factor, f , mean velocity, ue 5 ete.) 
depend on the following variables: d, , average flow depth; W, 
channel width; v , kinematic viscosity; k, » equivalent sand rough- 
ness; g , gravitational acceleration; S5 » bed slope. The depend- 


ence of E. on the other variables can be expressed as: 


E = A, ld 


: W,; 8, Ss EI (6.6) 


*? 


Without loss of generality, Equation 6.6 can be rewritten as: 


BE, = Ajldys 


y W, Vy ko» Uy,» g] (6.7) 


where u, = vgRS and R is the hydraulic radius. It is noted that 
for rectangular channels, the hydraulic radius, R , is defined if d, 


and W are specified. Dimensional analysis then gives: 


E. W k. ud) u, 
= A. j—, =, * (6.8) 
u,d 3 A d, Vv ore 
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It can also be shown that: 


k ud u 
* 
iy |e (6.9) 
* * ved, 


where f = 8(u,/U,)° and u,/Vgd, is a form of Froude Number. For 
fully rough turbulent flows, the friction Reynolds Number, u,d,/v 

has a minimal effect on E/uyd, and f, and u,/7gd, does not play a 
significant role in subcritical open channel flows and in ice-covered 


flows is not pertinent. Equation 6.8 and 6.9 reduce to: 


“5 W ae 
PEN OS ROE Ty sooo} 
k 
W s 
and: f = A. I>s = (6.11) 
6 E =| 


Combining Equation 6.10 and 6.11 also gives: 


- W (6.12) 
- Lalo mk 
u,d, A Es qj 


The normalised transverse diffusion coefficient, therefore, depends on 


the aspect ratio, W/d, , and the relative roughness, k/4, (or fric- 
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tion factor, f), for fully rough turbulent flows. 


If the relative roughness k 5/4, remains unchanged then 
E,/u,4, should approach a constant value for small aspect ratios. For 
very large W/d, » the width ceases to be an important parameter and 
should no longer influence the average size of the eddies. It follows, 
therefore, that if ko/d, is held constant, EO /uydy should approach 
a limiting value for large W/d, and for fully rough turbulent condi- 
tions. In subsequent ~-discussions, Equation 6.12 will be referred to 
instead of Equation 6.11 because the friction factor, f , is given 


for most of the available data. 


All the available laboratory data have been summarised in 
Table 2.1 and with the present experimental results, are shown plotted 
plotted in Figure 6.8. The plot reveals some apparent dependence of 
E/u,d, on W/d, as reported by Okoye (1970). However, it is evident 
that the normalised diffusion coefficient varies over a wide range at 
small W/d, and this is due to varying f. That this is true can be 
seen from the data of Miller and Richardson for which there was a 400% 
wotretion in the trictionstactor, f. Values of the normalised diffu- 
sion coefficient for W/d, > 30 are sparse and the available data for 
this ere are due to Okoye. Okoye's data for W/d, > 30 were obtained 
for smooth bed conditions only. Since f did not remain constant in 
all the reported data it will be erroneous to conclude that E/uyd, 
depends solely on the aspect ratio W/d,. It is believed that the 


variations observed in Figure 6.7 are primarily caused by varying f 
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and u,d,/v values. The hydraulic radius, R , is commonly used to 
describe channel geometry and thus it facilitates comparison of results 
from rectangular and non-rectangular prismatic channels. One could, 


therefore, replace d, with R in Equation 6.12 to obtain: 


A W 


u,R 12'R f1 Shing 


It must be remembered that Equation 6.13 is valid for fully rough tur- 
patent flows only. A plot of E,/u,R against W/R for all the avail- 
able data is presented in Figure 6.9. The effects of the friction 
factor, f£ , and the modified aspect ratio W/R cannot be separated 
since both parameters varied over a wide range in all the reported data. 
The only exception is the data of Miller and Richardson (1974). Their 
experiments were performed for a constant W/R and fully rough turbu- 
lent flow conditions. Their data can, therefore, be described by the 


relation: 
ete’ | haly FE ae (6.14) 


The plot of E,/u,R against f for the data of Miller and 
Richardson is given in Figure 6.10. It is seen that the normalised 
diffusion coefficient increases with an increase in the friction factor 


f when W/R is held constant and the flow is fully rough turbulent. 


The results of the ice-cover experiments do not fit the trend 
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of the other data when E/u,4, is plotted against W/d, (see Figure 
6.7) but they fit in well with the general trend of the open channel 
data when E/u,R is plotted against W/R. Leas tent fram Table 
6-3 that the normalised diffusion coefficient, E/u,R » for ice-cover 
conditions is, in general, only slightly less than the corresponding 


coefficient for open-water conditions. 


6.1.6 Attenuation of Peak Concentration 


Very close to the tracer source the tracer concentration at 
various flow levels will exhibit different decay rates along the axis 
of the tracer cloud depending on the position of the source in the 
vertical. For the case where the tracer is well mixed in the vertical 
it was shown in Section 3.1 that the decay rate of the depth-averaged 


peak concentration is given by: 


(3.15) 


where p is an exponent with a theoretical value of - 0.50; Cin is 
the depth-averaged peak concentration relative to the concentration for 
fully mixed conditions, and X, (= X/d,)) is dimensionless distance in 


the longitudinal direction. 


Figure 6.11 shows plots of Cin versus x for Runs 103, 


105, 201 and 205. The slope of each of the straight line plots is equal 


to the exponent p in Equation 3.15. A summary of the P values is 
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TABLE 6-4 SUMMARY OF DECAY EXPONENTS FOR THE ATTENUATION OF 
DEPTH-AVERAGED PEAK CONCENTRATION 


Decay 
Exponent 


Average -0.544 Average -0.542 
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given in Table 6-4. The results do not exhibit any definite trend 
regarding the dependence of p on any particular hydraulic parameter, 
because of the comparatively limited range of the parameters in these 


experiments. 


The plots, nevertheless, clearly show that for any two ex- 
periments performed at the same flow depth and discharge but with one 
having an ice-cover and the other a free water surface, e.g., Runs 
203 and 103, a longer distance is required to achieve the same value 


Ore eC This is in agreement with the earlier findings which estab- 


1m* 
lished a substantial reduction in E. due to the presence of an ice- 


cover. 


The average values of the decay exponent were found to be 
- 0.544 and - 0.542 respectively for open and ice-cover flows. 
These values indicate an attenuation rate which is slightly greater 
than the predicted value of p = - 0.50 for one-dimensional transverse 
diffusion. The reason for this is that the vertical distribution of 
tracer concentration along the channel axis is not strictly uniform 
until after an initial reach length which must depend on the flow depth. 
Within this reach the attenuation of maximum concentration is caused 
by Pere transverse and vertical diffusion. It can be shown, Czernuszenko 
(1973), that in a stream of large depth where the influence of the 
water-level and the boundaries on the tracer spreading are not marked 
the variation of C with X assuming Gaussian distributions in 


1m 1 


both the y and z directions is given by: 
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where Ey is the vertical diffusion coefficient. The diffusion in 

the y and z directions separately cause Cin to decay at a rate 
proportional to Soe SEY The fact that the average value obtained 
for p was close to the predicted value indicates that downstream of 


the first section of tracer measurement the effect of vertical dif- 


fusion near the axis of the tracer cloud was not very pronounced. 


6.2 Meandering Flume 


6.2.1 Concentration Measurements 


Transverse distributions of the normalised depth-averaged 
concentration, e/C, » for some test runs are presented in Figures 
On123a tof. Cs represents the tracer concentration when fully mixed 
with no tracer losses and is related to the tracer flux across any 


section by: 


N 
Ca -: Jouca: -[ 3 Deu, ¢, ta] [0-2 
j=1 
“| 


where c and u are depth-averaged concentration and velocity respec- 
tively; Q is total discharge; D is local depth; a is fraction of 


tracer which is lost and N = 26. 
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As could be expected, the data show that the tracer spreads 
much faster for centerline release of tracer than for side release. 
Figures 6.12a and f indicate a distinct displacement of the zone of 
peak depth-averaged Be Se eee the channel centerline. This 
displacement could be caused partly by the movement of the high velo- 
city filament from the inside to the outside bank, and partly by the 
spiral motions. It is observed, for example, to be much larger in 
run 401 than in run 404 and this seems reasonable because the strength 
of the spiral motion is greater in run 401 than in run 404 (see Figure 
5,26). A similar shift of the zone of maximum depth-averaged concen- 
tration can be seen in the transverse concentration distributions for 


run 302 (left-bank injection). 


The normalised concentration distributions in general appear 
to be smooth. For open channel tests patches of water having a high 
tracer concentration interspersed with patches having a low or even 
zero concentration were visible to the naked eye for some distance down- 
stream of the source. However, the fact that a reasonably smooth dis- 
tribution was always obtained at the first section of tracer measurement 
demonstrates that the sampling period used was long enough to smooth 


out most of the concentration fluctuations. 


Typical transverse concentration profiles in the first bend 
of the ice-covered channel at different levels and for various source 
locations are shown in Figure 6.13a. The plots indicate that for left- 


bank releases of tracer, the concentration distribution over a vertical 
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follows the pattern of transverse velocities in the vertical, the high 


concentration fluid being located near mid-depth in a vertical. 


The distributions for right bank and central releases of 
tracer are also consistent with what would be expected to result from 
the pattern of transverse velocities. This means that tracer is not 
uniformly mixed in the vertical as is usually assumed and that distances 
required to achieve vertical mixing are comparatively longer in the 
presence of transverse velocities. The observed lateral displacement 
of the tracer cloud in run 401 is attributable to the fact that the 
tracer was released at oT = 1/6 for this particular run and it is 
evident from Figure A9a that significant transverse velocities exist 


at this section. These velocities could initially shift the entire — 


tracer cloud laterally. 


Corresponding transverse concentration distributions at var- 
ious levels and for different source locations are presented in Figure 
6.13b for some open channel tests. Again the concentration distribu- 
tions are consistent with the measured transverse velocity distributions. 
Run 302 indicated rather large differences in concentration over a 
vertical and a marked lateral displacement of the zone of high concen- 
tration. It is noted that the average velocity for this run is 0.80 
ft/sec and the maximum observed transverse velocity near the water sur- 
face between 8, = 0 and T/2 is about 0.13 ft/sec. If a vertical 
line source is imagined to be located at the left bank then in the 


absence of turbulence the observed maximum w would cause a lateral 
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TRANSVERSE CONCENTRATION PROFILES AT VARIOUS FLOW 
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FIGURE 6.13b TRANSVERSE CONCENTRATION PROFILES AT VARIOUS FLOW 
LEVELS AND AT SECTION 95 = 7/2 
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spread of tracer of about 2.25 ft (or z/w = 0.9) at 05 Ties 161s 

lateral spread appears to be comparable to the observed spread at y/D 
= 0.75. Again in the absence of turbulence there will be no lateral 

spread in the lower half of the flow due to the effects of transverse 
velocities. The effect of turbulence is to even out any vertical dif- 
ferences in concentration, however, the large differences observed in 
run 302 are probably due to the fact that vertical mixing takes place 


over a relatively large depth. The differences for runs 304 and 306, 


Figure 6.11b, are seen to be less because of the smaller depth of flow. 


Turbulence causes the observed vertical differences in con- 
centration to diminish with distance downstream. However, if these 
differences are still large when the flow enters the second bend the 
subsequent reversal of the spiral motions will result in a gradual 


modification of the vertical concentration profiles. 


The root mean square concentration deviation, oo » often 
termed the coefficient of variation, of the depth-averaged transverse 


concentration distribution can be defined as: 


1/2 
[e, - c,}7 


C_ is a good indicator of the degree of transverse mixing and is anal- 
ogous to the non-uniformity coefficient, 6 , used in Chapter 5. 


Figures 6.14a and b give the longitudinal variation of Cy for all test 
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Coefficient of Variation (Cy) 
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FIGURE 6.14a COEFFICIENT OF VARIATION CURVES FOR ICE-COVER TESTS 
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FIGURE 6.14b COEFFICIENT OF VARIATION CURVES FOR OPEN CHANNEL 
TESTS 
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runs. 


It is to be noted that low values of C. imply a high degree 
of mixing. As expected, much better mixing was always obtained for 
tracer introduced on the channel centerline compared with the side re- 
lease of tracer for the same flow conditions. The diagrams also reveal 
a tendency for the rate of mixing to have an anomalous decrease at the 
beginning of the second bend (8, = 0 to oD) = 1/3). This is clearly 
demonstrated in runs 405 and 408, the former even indicating an increase 
in Le This phenomenon is presumably associated with the decay and 


reversal of the spiral motion from the first bend. 


There is evidence from Figures 6.14a and b that the tracer 
undergoes different rates of mixing within the same bend depending on 
the the source position. This is consistent with previous observations 
in Chapter 5 that the spiral motion varied both laterally and longitu- 


dinally. 


It is appropriate here to compare some of the above results 
with that of a straight flume. If it is assumed that BE /uydy = 0.20 
for straight open channels, then according to the procedure of Ward 
(1973) the distance, Xn , required to achieve a given percentage degree 
of mixing (defined as [1 - Cy] x 100%) is given by xX /W s T,/0.20 ° 
Uta W/d, Ty is a variable which depends on the source position 
and the specified degree of mixing. Run 304, for example, achieved 


approximately 95% mixing at 0, m/f 3 According to Ward's method an 


equivalent straight channel with the same W, d,, U, and u, as in run 
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304 would require at least 6 times this distance to achieve the same 


degree of mixing. 


6.2.2 Evaluation of Transverse Exchange Coefficient 


6.2.2.1 Method of Moments: From the theoretical considera- 


tions of Chapter 3, the following equation was derived for evaluation 
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To assess the importance of the dispersion term u'c' a 
straight line was fitted by least squares to the concentration data in 


a vertical and the c' 


values were determined although the vertical 
concentration profiles were sometimes far from linear. The velocity 
distributions over verticals with concentration data were determined 


by interpolation and the u' values computed. The u'c' 


values were 
then evaluated by numerical integration. It was found that u'c' was 


in general about two orders of magnitude smaller than uc and, there- 


fore could be neglected. 


By a suitable choice of coordinate system it is generally 
possible to make the depth-averaged transverse HRS, w , negligi- 
bly small and for this situation the convective term wc becomes 
unimportant. Such a coordinate system was employed by Yotsukura and 
Cobb (1972) who used the cumulative discharge, q, » as an alternative 
transverse coordinate. This system of coordinates will be utilised in 
Chapter 7 to analyse the field data and it is shown there that this 
method gives practically the same value of ee as that obtained by 
using the meandering coordinate system. This indicates that we is 
of secondary importance compared with wile In the laboratory, the 
magnitude of w was of the same order as the error in measurement of 
w hence wc could not be accurately determined. However, because 


of its apparent unimportance in the more extreme field situation it was 


assumed so in the laboratory. 


The major effects of curvature on the mixing, therefore, must 
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be caused by the dispersion term w'c'. A general comparison of the 
relative importance of w'c' and transverse diffusion can be made 
using, for example, the data from run 401 (centerline injection). It 
is noted that the average velocity for this run is 0.85 ft/sec and the 
average relative transverse velocity between tracer particles located 
in the central and either the upper or lower layers of the flow is 


about 0.085 ft/sec. The transverse displacement rate due to this re- 


0.085 
0.85 


length. An estimate of the diffusion rate can be obtained by assuming 


lative transverse velocity would be 


= 0.1 ft per foot of channel 


that the concentration distribution due to diffusion alone is Gaussian 
and that BE /u,4, a 8 Pat 8 TE it is noted that in’ run 401 Uo/ ts = 


12.65, the above assumptions indicate that the standard deviation: 


2E, u a2 


would grow at a rate of about 0.029 ft per foot length of channel. The 
displacement rate due to transverse velocities is, therefore, much 
greater than that caused by diffusion. The same was observed to be 
true for the open channel tests. It is clear from the above illustra- 
tion that the transverse spreading of tracer could be reduced apprecia- 


bly in regions where the pattern of transverse velocities reverses. 
The function @$ generally varies in an unknown manner. In 
the present analysis it is assumed equal to unity. For this situation 


ES = E. and since u'c' << uc. Equation 3.27 can be simplified to: 


dM, (X) 
dX 


- G(X) EO ° B(X) (6.16) 
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in which: 
Wo Wo 
M(x) = | Ducn dn | Ducdn (6.17) 
yal 7 


and represents the second moments of the tracer flux about point of 


tracer release: 


W 
: es A ee Ni: 
2 | Das - —)(we +t w'e')dn 
c 
Le 
(ey = SSS (6.18) 
2 
Joucan 
oT 
and: 
W W 
2 2n + Z 2 nN + Z5 \D 
| peu -£ ydn + | end = > ) an 20 
W . Ws 
B(X) =2 = 
2 
[oucan 
ud) 
W, . 
9_ tpen( + °)}4 
| an {Den(1 + z )tdn 
7 
ee) (6.19) 
Woe. 
foucan 
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(21.8) 


ir Pe ae shi ¥ 


30 | ; 


Ke" Bes; ef ot y 


Vga 
as a 


bu neal 
(21.3) | 


171. 


Equation 6.16 simply states that the change in second moments of the 
tracer flux results from the spreading of tracer caused by transverse 
velocities and diffusion. The transverse exchange (or diffusion) 
coefficient E. can be determined using Equation 6.16 provided the 


variations of M, (X) , G(X) and B(X) are known. 


In practical situations we are more often interested in the 
overall mixing due to both turbulence and transverse velocities and 
it is convenient to define an exchange coefficient which accounts for 
these two factors. It was shown in Chapter 3 that the relevant equa- 
tion for the above situation is Equation 3.29. Again, if as is 
regarded as equal to unity and therefore ae = ee then since 


uc >> u'c' Equation 3.29 can be simplified to: 


E ow * B(X) (6.20) 


The variation of the second moments of the tracer flux, M, (X), 
with longitudinal distance, X , are presented in Figures 6.15a and b, 
for some test runs. In general M, (X) varied non-linearly with X 
and this non-linearity is believed to be related to the growth, decay 
and reversal of the spiral motions. Some of the plots indicate a 


flattening out and even a decrease of M, (X) between 8, = 0 and 


n/3. 
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There is a significant dip in M, (X) for run 406 at 8, = 1/2. 
A study of velocity profiles for this run in Figures 5.13 | 
reveals that the warping of the plywood (ice-cover) caused the maximum 
flow depth to occur at the inside of the bend instead of the outside 
as could be anticipated. The change in cross-sectional shape between 
8, = 1/3 and mt/2 could cause a net transverse flow towards the right 
bank and thus inhibit the spreading of the tracer to the left bank. 
This will ae oct itself by a reduction in the observed value of 
M, (X) and probably by a slight increase in Cy (see Figure 6.14a). 
If the above physical reasoning is correct then there must be a cor- 
responding increase in M, (X) at 8, = 1/2 for left-bank release of 


tracer under the same flow conditions. This is confirmed by the varia- 


tion of M, (X) for run 405 in Figure 6.15b. 


As the tracer concentration becomes more uniform across the 
channel, the second moments of the tracer flux, M, &) » should approach 
a constant value. This is illustrated by runs 304 and 401 in Figure 
6.15a, b. From Equation 6.19, the following functions can be defined 


for o, = 1.0: 


B, (X) = 2 ora + ———*) dn Jo ucdn (6.21) 
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Wy W, 
= ye 7 oad Se) 
B, (X) = 2 jena te ares mn dn [p fine dy (6.22) 
¢ 
Wy Wy 
L 
W W 
i 2 
5 € n+ z) 7 for. 
B, (X) = 2 ES {Den (1 aaeaaee [p UNGHGan) (6.23) 
W e W 
1 1 
so that: 
B(X) = B, OD) + B, (X) - B(x) (6.24) 
The longitudinal variation of Bi» Bo» B3 and B are given 


in Figures 6.16a, b and c for some test runs. It is evident from these 
plots that for centerline release BL remained reasonably constant 
over the length of the test reach. However, for injection at the right 
bank, B, was initially high but decreased to a fairly steady value 


- 


downstream. When tracer is introduced near the left bank By increases 


/ 
to a maximum before dropping to a near constant value. These differences 


in the observed By distributions for the various source locations can 
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be explained by imagining a vertical line source of unit width and 
uniform concentration oe With its centroid located at_,z . itvis 
o 


to be noted that tO = 9 ft and that within the first bend: 


W 
: Ss 2n Z 5 
By = | De(1 + — % 
il 


where Q, is the tracer flux. The following relations can be obtained 


for 5 at the section of tracer release (X = 0) for the various source 


i 


locations: 


Source at left bank, 4 = +1.20 B, (0) = 0.867 DC ./Q, 
Source at centerline, Dh 0 B, (0) = 1.0 DC./Q, 
Source at right bank, Z3-4 1.20 B, (0) = 1.133 DC ./Q, 


The above relations predict an initially high B, value for injection 


il 


at the right bank and a low value at the left bank. These differences 
in By values therefore appear to be largely due to the system of 


coordinates used here. 


B, expresses the importance of depth variations within a 
cross-section on the mixing. Such variations in the laboratory experi- 
ments result from superelevation and unevenness of flume bed and ice- 


cover. It is evident from Figure 6.16 that variations in depth were 


generally insignificant in the laboratory tests. 


B3 does not influence or affect the computed value of E.. 


until a significant portion of tracer encounters both banks. For the 
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situation where the source position is near the side, B, is expected 
to decrease slightly with distance downstream because of the decay of 
concentration along the bank on which the source is located. However, 


as more tracer encounters the opposite bank B, might increase and 


3 
approach a constant value which is attained when the tracer becomes 
uniformly mixed laterally (see Figure 6.16c and runs 302 and 306). 
Strong variations in By were observed for right bank injection com- 
pared with left bank tracer release. Figures 6.12b and e, for Sehneae 
reveal that the concentrations at the right bank are much higher when 
tracer is released near this bank compared with the corresponding case 
of left bank injection. For source located at the channel center B, 


will be zero until the tracer reaches the sides and then will increase 


until tracer is uniformly mixed. 


From the above discussions B (= Bi + B, - B,) is expected 
to decrease gradually with distance downstream of the source and to 
approach zero as the tracer becomes fully mixed. This trend can be 
seen in Figures 6.16 for some of the test runs. Some test runs, for 
example 403 and 406, do not exhibit this expected variation in B and 
this is attributed to the fact that these runs did not achieve any 


appreciable degree of mixing even at o, = 21/3. 


The longitudinal variation of B(X) is presented in Figures 


181 


6.15a and b together with that of M, (X) for some test runs. Numerical 


values of M, (X), B(X) and 2. for all test runs are summarised in 


Table 6.5. To facilitate comparison of ice-cover experiments (series 
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400) with those having free water-surface (series 300), tests performed 
under similar flow conditions are given the same number in their respec- 
tive series. Comparison of the Bie values for the open channel and 


ice-cover tests will be given later. 


6.2.2.2 Integral Method: The integral method is used not 
only to evaluate the average transverse mixing coefficient over a 
channel reach but also the variation of E ow over transverse portions 
of the same reach. It is to be noted that E ow in this case repre- 
sents the average transverse mixing coefficient for either the whole 
or part of the reach. From the theoretical considerations in Chapter 


3 it was shown that: 


z Zz 
ad. Bees A TE = 
| a | D(u c + u'c')dz]{ dz = 
Wie 
Jy. Vso ) 
E Dhic ~ c —~—— dz (3.36) 
1 Oz 
ok een 


The integrals in the above equation were determined numerically by per- 
forming the integration from the bank of tracer release for progres- 
sively increasing or decreasing values of z until the entire width 

of the tracer cloud was covered. The value of the term in square 
brackets on the RHS of Equation 3.36 was evaluated by averaging the 


corresponding values at the beginning and end of a reach. The cross- 
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channel variation of Be within a reach is shown plotted at the down- 
stream section of the same reach in Figures 6.17a to d for some test 

runs. It is mentioned that the observed lateral and longitudinal varia- 
tions in ae are actually due to variations in the depth-averaged 
transverse mixing coefficient ee since Ew is considered to represent 


the average of the generally variable CLE 


As the integrations in Equation 3.36 proceed from one bank 
to the other the effect of small variations in [ee on Ew are not 
discernible from the plots and appreciable changes in ao only cause 


a slight change in Eves especially when the integration covers a large 


percentage of the channel width. 


Figure 6.17 shows that negative values of pe occur near 
the left bank when the source is located close to this bank. This can 
be seen in both ice-cover and open channel tests and for some distance 
downstream of the source. This feature is not unreasonable because the 
spiral motions could initially displace the zone of maximum depth- 
averaged concentration away from the left bank (see Figures 6.12d and 
e) and at the same time transport low depth-averaged concentration fluid 
from the sides against the concentration gradient to regions of high 
depth-avetazed concentration consequently causing a negative ae As 
expected, this does not occur for tracer released near the right bank. 
The plots also indicate low and even negative ps values at the be- 


ginning of the second bend and this accords with the results from the 


moments method. 


ee! rn | je pe . 


sy = 


~aeod edt 28 beodett guetta ek dTasot maaan A 20 : 
ses? enor xo? b of sees eprugtt chennai 

-stxev Lsekbugignol bak axeost bovrando ‘arta ‘a lair 

begazeve-dageh 9d3 at.gnohipt sav oF sub eLinae 

inave1get 03 batebtemon at a sonke a; teakottieos — 

“al? stay Sie ody 0% Ave 

poe Sa of : 


ined ado ec 22 bosvong ‘see potsaupt nt ‘yaotaesgeamt ods ae 
y hy eM 


a wat: “eb aaotantaae CTann Yo 2o0te a eae 


, 


saueo tine ee nt sognad> sidetoouage bra sate cea idiaiea 
ui) pe 


sgzel 6 exavas: ‘natsorgeoal fe amie hinting pet anes: 


en ob 


rao me Se esutey ounseiied ‘gues aworla XEa , ; 
mao ahtt ,slend absle om satis heteoal’ ef saxo0e oats, 2 
ostaselh anne Tol bie esesd beers 3 staal 
eth sedond sldeworssuay: toa iim axygeot abst sores as heh: “a " 
- tiget ilmbiena bo pret vine aoalqesd vile i be ae snot som 
baw i61.9 sowgit soe) dats Aetna mont (win sotgnzansonos bs | 
btu notssrsagsn03 cations snoganess « ents ius - = oan 
seit Ye saeary ot. 3 berg ha) 2660: 
an “ane avigegou. © whale $ 
"ci aia <c eees Se ie 
~oé oi de Seuley _entsnyaiobva te Wot Banatbad onte wake at 
sila mnt ‘ gi Tueas 7 v — elds boa alse 36 satinia 


, : ' i. y ; - 2% * 
: (a ed) a ne kan ners ees 
*ie F. —_ - 
3 3 a lhe r : i : > ! Me a x 7 —@ “A 
- a on $ tal j 
i | se J . 
» 7 -§ + i 7 1 “5 
| i i i. qu 
_ ‘ ae | é 7 ae AG 


Run 405 (LB) 


Exchange Coefficient, Ezw (ft /sec) 


z/W 
(a) 


Run 403 (RB) 


184 


FIGURE 6.17 LONGITUDINAL AND LATERAL VARIATION OS THE MIXING COEFFICIENT 
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For situations where source is located at the right bank Ew? 
or more correctly ee » is again seen to exhibit marked variations near 
this bank in both bends. Figures 6.17 show that there are similarities 
in the general pattern of variation of EAs for the same source posi- 
tion but different flow conditions. This is dramatically illustrated 
by runs 306 and 403. Both runs indicated a high hee value near the 
right bank between 0, = 7/6 and mt/3 although the pattern of spiral 
motions was different in both runs. No explanation can be offered for 
this behaviour in EB, but it is speculated that it might be due to 
some strong disturbance within the flow in the neighbourhood of the 


right bank between 8, = 17/6 and m/3. 


It is evident from Figures 6.17 that the depth-averaged trans- 
verse mixing coefficient varied both along and across the channel. The 
variation of aes for the various test runs is listed in Table 6.5 


together with the values obtained by the method of moments. 


The integral method was found not to work for situations 
where the source was located at the channel centerline. It yielded 
values of E ow which were on the average an order of magnitude greater 


than those given by the method of moments. 


It is evident from Table 6.5 that for similar conditions of 
discharge and flow depth and for the same source position the mixing 
capacity as measured by a5, was of the order 3 to 6 times greater for 


open channel flows than for ice-covered flows. The straight flume 
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experiments indicated that the diffusion coefficient, E. » for open 
channel flows was in general twice as large as that observed for ice- 
covered flows for correspondingly similar flow conditions. It is not 
unreasonable to assume that this observation is also true for the 
flows in the meandering flume because the flume bed and ice-cover 
roughnesses were the same as in the straight flume tests. Since the 


mixing coefficient, es » represents the rate of mixing due to both 


turbulence and bend-generated transverse velocities the further increase 


in pee (1.0 to 4 times, relative to Be for ice-covered channel 
flows) for the open channel tests can be attributed to additional mix- 
ing caused by transverse velocities. That this is true can be seen by 
comparing Figures 5.6 and 5.10. The average values of the strength of 
spiral motion, S , for open water tests are observed to be about 1.5 
to 2.0 times greater than the corresponding ice-cover tests. The 
comparatively high values of ‘geod for open channel tests, therefore, 
seem to be reasonable. One factor which must be noted is that the 


curved flow could modify the structure of the turbulence field, but 


there is no information available on this modification. 


6.2.3 The Normalised Transverse Mixing Coefficient 


A summary of the normalised transverse mixing coefficient 
k (= Eg! Ux) for all test runs is given in Table 6-6. The values of 
k were computed using the overall average hydraulic radius, R , and 


the average shear velocity, u, , and are shown plotted against X in 
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Figures 6.18a and b. There is a tendency for the Integral method to 
give higher values of k in the first bend especially at the first 
section downstream of the source and lower k-values in the second con- 
pared with the moments method, when tracer is released from the left 
‘bank. The reverse tendency is seen to be true for the case where 
source is located on the right bank. In general, the values of k 
obtained by both methods are remarkably close and both methods exhibit 


similar patterns in the longitudinal variation of k. 


The discussions which follow are based on results obtained 
by the method of moments but the same general comments also apply to 
results obtained by the integral method except for situations where 


tracer was released from the channel centerline. 


It is evident from Figures 6.18 that the longitudinal varia- 
tion of k is influenced by the source position and probably reflects 
the role of transverse velocities in the mixing. As could be expected 
the normalised transverse mixing coefficient for centerline release of 


tracer was observed to be greater within the first half of the first 


bend. k reaches maximum values for different source positions between 


95 = 7/2 and 7 of the first bend while minimum k-values occurred at 
the beginning of the second bend. It is of interest to note that the 
decay and reversal of the spiral motions also occurred, in general, 


between 9, = 0 and m/3. 


High k-values were observed in the second bend when tracer 
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FIGURE 6.18a LONGITUDINAL VARIATION OF THE NORMALISED TRANSVERSE 
EXCHANGE (MIXING) COEFFICIENT 
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was introduced at the right bank. This is probably due to the fact 
that when the source is located at this bank the spiral motions in the 
first bend appear to inhibit spreading of tracer in that bend while 

in the second bend they seem, initially, to enhance the spreading of 
tracer. Evidence for this is seen, for example, in Rie res: Gals and 
runs 403 and 406. Both runs indicated z sharp increase in second 


moments between oP = 7/6 and m/3. 


Negative values of the normalised transverse mixing coeffi- 
cient occurred for some test runs at the beginning of the second bend. 
This is not unexpected. Although this violates the concept of gradient- 
type mixing it is realised that the gradient-type mixing concept is an 
approximation and cannot be expected to apply to all flow situations. 

It is possible to obtain negative k-values in regions of curvature 
reversal as this is initially accompanied by a net lateral discharge 
towards the outside bank. This net flow can inhibit spreading of 

tracer released from the outside bank and produce a negative k-value. 

A similar situation applies when longitudinal changes in cross-sectional 
shape cause significant transverse velocities, for example, run 406 and 


between 0. = 1/3 and 1/2 ( Figures 5.13 and Figure 6.18b ). 


Further downstream in the second bend the normalised trans- 
verse mixing coefficient is again seen to increase to a maximum in 
the middle third of this bend before decreasing again towards the bend 
exit. The longitudinal variation of k for some ice-cover tests was 


rather erratic and is probably caused by the warping of the plywood 
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(ice-cover). This was more pronounced in runs 404/5/6. It is noted 
that the flow depth for these runs was comparatively small and conse- 


quently the results were more affected by the warping of the ice-cover. 


The tracer was observed to be almost fully mixed at o, = 7/6 
for run 302 which indicated high values of k in the first bend. The 
k-values evaluated for subsequent reaches beyond 8, = 7/6 for this 
run were found to be unreasonably large and this could only be traced 
to measurement errors. As discussed previously when tracer becomes 
more uniformly distributed laterally the variable B(X) tends to zero. 


Any errors in the measurements of velocity and concentration, therefore, 


can result in very unrealistic values of k. 


The values of k averaged for different source positions but 
for the same flow conditions are shown plotted against X in Figure 
6.19. It is evident from this figure that maximum values of the norma- 
lised transverse mixing coefficient occur in the middle third of both 
bends and minimum values occur at the beginning (0, = 0 and m/3) of 
the second bend. When the k variation in Figure 6.19 is compared with 
the variation of spiral motion strength in Figures 5.6 and ot ae ee 2 
evident that k seems to be strongly related to the growth-decay- 
reversal cycle of the spiral motion. With only two bends it was not 


possible to relate the probable cyclic variations of k to that of S. 


Table 6-7 shows a comparison between the normalised mixing 


coefficients obtained over the same reach for the same source position 
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Runs 401, 402, 403 


Runs 404,405,406 


Runs 407, 408 


Bend Angle (radians) 


FIGURE 6.19 LONGITUDINAL VARIATION OF THE AVERAGED NORMALISED 
TRANSVERSE MIXING COEFFICIENT FOR DIFFERENT FLOW 
CONDITIONS 
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TABLE 6-7 COMPARISON OF E og! tak FOR THE 
SAME REACH AND SOURCE POSITION 
AND SIMILAR FLOW CONDITIONS 


Ice-Cover Open Channel 
Tests Tests 


E 
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and similar flow condition in the open and ice-covered channel tests. 
The magnitude of k is seen to be considerably greater for open-water 
than for ice-covered conditions. Since an ice-cover also reduces the 
product u,R by a factor of approximately 2/2 » if depth remains 
constant, the results again indicate that an ice-cover can reduce the 


mixing capacity in channel flows substantially. 


6.2.4 Parameters Controlling E /u,R 
zw 


——— OW 


In turbulent meandering channel flows, the dependence of bulk 
flow properties such as the mixing coefficient, ES » average velocity, 


im » and friction factor, f , on other parameters may be expressed as: 


Ess, = 7, (QR, Wes ko: g5 S5° to L) . (6.25) 


where the variables have their usual meaning and L represents the 


Overall bend length (measured along the channel centerline). 


Dimensional analysis can be used to give: 


Ew IRE RF pd Ks og oe re Le (6.26) 
> 9 > > > > : 
u,R Zaha, Vv BR R ro 


where u,/vgR is a form of Froude Number. 
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Again for fully rough turbulent flows the friction Reynolds number 
u,R/v has a minimal effect on E og! Bak and the Froude number, u,/vgR, 
does not play a significant role in subcritical open channel flows and 


in ice-covered flows is not pertinent. For these conditions Equation 


6.26 simplifies to: 


E k r 
A jes Oh W s (id vo 
c 
or: 
E r 
Zw W CUPL 
u,R ‘ uP R* = rR Cr (6.28) 


The present data cannot be Meh to explore the importance of all the 
various parameters in Equation 6.27 or 6.28 because of the limited 
range of these parameters in the experiments. Nevertheless, the ice- 
cover data can be used to gain some insight into the dependence of 
Ey! Ua on the other four parameters. The values of the parameters 


in Equation 6.28 are summarised in Table 6-8. 


TABLE 6-8 NUMERICAL VALUES OF VARIABLES IN EQUATION 6.27 
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R Tr 
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90 | 5.23 | Rough 
404/5/6 | 0.658 | 36.8 | 0.0489} 132 | 5.23 |] Rough 
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It is evident that there is a tendency for k to decrease 
with increasing r/R and W/R. The effect of roughness is not too 


clear because the data is insufficient. 


Fischer (1969) investigated the effect of bend generated 
spiral motions on transverse mixing. He obtained an expression for 
the additional mixing caused by transverse velocities, that is w'c' 5 
based on Rozovskii's transverse velocity distribution in the vertical 
for fully developed turbulent open channel bend flow. By further 
assuming that w'c' = - Ey dc/8z , where e€, is the transverse dis- 


d 


persion coefficient, Fischer obtained: 


a Copel eo: =) = (6.29) 
ads u, ce 2 


in which kK is von Karman's constant and values of the function I 


are given by Fischer. Fischer recognized the limitations of the above 
equation and suggested that it should be used for order of magnitude 


estimates only. Equation 6.29 can be rewritten as: 


€ Wy AW 2 
oa ay kes x] (6.30) 
u,R u, To 


According to Fischer's predictions, therefore, €,/u,R should decrease 
with roughness. This appears to be consistent with earlier observations 
that the strength of the spiral motion decreases with roughness for 


similar conditions of depth and discharge. Also €,/u,R should decrease 
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with increasing r/R ratio. If it is assumed that r/R does not 
affect the diffusion coefficient there is some evidence in Table 6-8 

to support this prediction. It is noted that because of the assump- 
tions inherent in the derivation of Equation 6.29 the aspect ratio, 
W/d,. » does not appear in this equation. The importance of the various 


parameters in Equation 6.28 will be discussed further in Chapter 7. 
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CHAPTER 7 


FIELD INVESTIGATIONS 


Velocity and tracer measurements were made in the Lesser 
Slave River (Alberta) in August 1972 and February 1973, and these mea- 


surements are described and discussed in this chapter. 


7.1 Description of the Test Reach 


The field tests were conducted on the Lesser Slave River 
(Alberta), approximately eight miles downstream of Lesser Slave Lake. 
This site was selected because discharge from this large lake remains 
constant over prolonged periods, the water quality is suitable for dye 
tracing, reasonable depths of flow can be obtained under ice-covered 
conditions, and there are no tributaries within the study reach of 
12,000 ft. length. During the winter test the river width varied be- 
tween 90 and 160 ft and the average depth was between 5.5 ft and 8.0 ft. 
The maximum depth was found to be less than 14 ft during both winter 
and summer tests. The river has an irregular, but marked meander pat- 
tern with a sinuosity of about 1.8. The bed material is predominantly 
fine sand with a median size of 0.20 mn. The hydrologic and geomorphic 
characteristics of a nearby, similar reach of the Lesser Slave River 
are listed by Kellerhals, Neill, and Bray (1972). The Water Survey of 
Canada maintains a gauging station approximately six miles downstream 


of the study reach. A plan of the study reach is shown in Figure 7.1. 
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Bok Experimental Procedure 


Tego Velocity Measurements 


For the velocity and tracer measurements under ice-cover, 
holes were drilled through the ice at the six numbered cross-sections 
shown on Figure 7.1. Velocity measurements were taken at 1 ft inter- 
vals in the vertical and at 20 ft intervals across the river. A few 
more detailed vertical velocity profiles (0.5 ft intervals) were also 
obtained. Under open water conditions, velocities were only measured 
at 0.2 and 0.8 of the local depth, D , along verticals spaced at 10 


ft intervals at Sections A, C, and F shown in Figure 7.1. 
7.2.2 Tracer Measurements 


Rhodamine WT dye diluted with stream water to about 2% by 
weight was used as tracer. For the winter tests it was injected at a 
constant rate of 13.5 eee just below the underside of the ice at a 
point approximately 60 ft from the left-bank of the river at Section 
0. About 90 minutes elapsed after commencement of injection until a 
steady state concentration distribution was established at Section l. 
After this time, sampling was begun at Section 1 and progressed gradu- 
ally downstream to Section 5 where it was completed 6 hours after the 


start of the experiment. 


A sampler, capable of taking 60 cc samples simultaneously at 


various depths up to 8 ft,.was used. Samples were collected at 3 ft, 
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5 ft and 7 ft below the free water surface, which corresponds to depths 
of approximately 1.5 ft, 3.5 ft and 5.5 ft below the underside of the 


ice, and at 10 ft intervals across the channel for each section. 


For the open-water test, the tracer was injected 15 ft from 
the left bank at Section A, and samples were taken at Sections B, C, D, 
and E at depths of 1 ft and 4 ft using a 10 ft spacing across the 


channel. 


7.3 Experimental Results 


i20.L. Lee, Cover 


Ice-thickness was found to vary between 1.0 and 3 ft, being 
greatest near the banks except where the zone of maximum velocity was 
close to the bank (Figure 7.3). Three slabs were cut from the ice 
cover and turned over to reveal the relief on the ice underside (Figure 
7.2). A surprisingly regular dune-like pattern is evident in Figure 
7.2(a). The ridges are perpendicular to the main flow direction. 
Figure 7.2(b) shows a complex, wavy and irregular pattern, while Fig- 
ure 7.2(c) shows a relatively smooth underside. With only three 
samples, no firm conclusions concerning the spatial HTS eIbietontot 
the 3 patterns can be drawn but a general smoothing of the ice under- 


side in the downstream direction was indicated. 


7.3.2 Hydraulic Data 


The hydraulic data for the winter test are summarized in 
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FIGURE 2. RELIEF ON ICE UNDERSIDE 
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Table 7-1 and velocity contours for the six cross-sections are presented 
in Figure 7.3. The point of maximum velocity was observed to be located 
anywhere between 0.2D and 0.7D, measured from below the underside of the 
ice cover. However, it generally occurred within the upper half of the 
flow. The vertical location of this point of maximum velocity is con- 
trolled by two main factors: (i) the relative roughnesses of ice- 
cover and river-bed and (ii) the spiral motions induced te the bends. 
Some of the velocity contours of Figure 7.3 indicate the presence of 
secondary current cells but the exact number and extent of the cells is 
generally not clear. The variation in discharge over the three days on 
which velocities were measured is probably caused by differing wind 
conditions in the region of open water at the outlet of Lesser Slave 


Lake. 


The hydraulic data for the corresponding summer test are 
listed in Table 7-2 and two observed lateral distributions of mean velo- 
city (for a vertical) are shown on Figure 7.4. Exact correspondence 
between summer and winter cannot be achieved because winter flows are 
generally lower than summer flows, and at equal flow, winter stages are 


considerably higher. 


The slope of the free water surface, Sy » was observed to be 
0,095 <= 10> under winter conditions, and this should correspond close- 
ly to the slope of the energy line and to the open-water slope. Lt is 
close to the value of 0.11 x ign listed in Kellerhals, Neill and Bray 


(1972) for the nearby reach of river. 
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The average values of R and u, are 3.22 ft and 0.10 ft 


x 
a an 
s for ice-covered conditions and 8.3 ft and 0.16 ft 5 z for open- 


water conditions. 
7.3.3 Tracer Distributions 


The observed transverse tracer distributions are shown in 
Figures 7.5a and b for winter and summer respectively. It is to be 
noted that Cy is the concentration for fully mixed conditions with 
no tracer losses. The effects of secondary currents due to bends are 
noticeable in both winter and summer data. The winter data in parti- 
cular show the zone of highest concentration moving downwards (between 
Sections 1 and 3) and then upwards (Sections 3 to 5). The displacement 
of the zone of maximum concentration with the shifting thalweg can be 
observed between Sections 2 and 4. The distances needed for vertical 
and transverse mixing under ice-cover conditions appear to be of the 
same order of magnitude in the Lesser Slave River. This differs from 
the case of reasonably straight broad channels, where the distance 
for vertical mixing is much shorter than the distance for transverse 


mixing. 


ree Computation of the Transverse Mixing Coefficient, E = 


No direct field measurements of the transverse velocities, 
w , were made in this investigation and their effects must, therefore, 


be included in a modified exchange coefficient, ae » defined in 
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FIGURE 7.5b TRACER DISTRIBUTIONS, SUMMER. 
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The relevant equation for evaluating Ew is Equation 3.30 


and it is to be remembered that this equation is valid only for rivers 


with zero depth at the banks. Again, 3 is assumed equal to 1.0 so 


that Equation 3.30 simplifies to: 
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was fitted by least squares to the concentration measurements in a ver- 
tical and the c" values were determined for the points where u' was 
known. The velocity distribution over verticals on which no velocities 


were measured were found from the velocity contours. The values of 


"c' exceeded 3% of the corresponding uc value on only 3 verticals 


u 
out of 65 and these were located near the tails of the transverse 

tracer distributions where the uc values are themselves small. For 
practical purposes u'c' can, therefore, be dropped from Equation 755. 
This does not mean that the vertical concentration gradients are unin- 


portant, since their major effect would appear in the term w'c' which 


has been included in E ow’ Equation 7.5 now reduces to: 
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= —E B(X) Cie7) 


The definitions of M, (X) and B(X) can be obtained by comparing 
Equations 7.6 and 7.7. Since the width of the river varied over a wide 
range, the calculated second moments had to be adjusted to allow for 


this variation. The formula used was (Fischer, 1967): 


i Pi 
in which (My), is the corrected second moment; (M) . is the actual 


or calculated second moment; Ws is the maximum width at a cross- 
section and ae the average maximum width for the test reach. The 
adjustment of M, (X) as suggested by Fischer is not entirely unrea- 
sonable because convergence or divergence of the flow does not neces- 
sarily imply reduced or increased mixing. The values of (Mo) o> (Mo) o> 
B ata Ew are tabulated in Table 7-3 and in Figure 7.6, plots of 


(M, ) and B against X are presented. It is evident that the 
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decrease in B with longitudinal distance, X , follows the same 
trend as was observed for the laboratory experiments which were con- 


ducted at approximately constant depth. Possible lateral variations 


of the two terms: 
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TABLE 7.3 NUMERICAL VALUES OF SOME VARIABLES 
IN EQUATIONS 7.7 and 7.8 
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are illustrated in Figure 7.7. It is evident that depth variations 
modify the distributions of Fi and Fy but when Fy and F, are 
integrated over the cross-section, the effects of depth variations on 

B and consequently on the normalised mixing coefficient, k , are less 
obvious. Nevertheless, some understanding of the importance of depth 
variations on the mixing may be gained by comparing the results from 


channels of non-uniform cross-section with those from equivalent rec- 


tangular channels. This aspect of the mixing is discussed later. 


The average value for the transverse mixing coefficient, 
E ow » over the length of the test reach is 0.113 fea: for the winter 
test and the normalised mixing coefficient, k , is 0.36. The open 
water values are Ee = 0.433 Fea and k= 0.33. Although this 
k-value is only slightly less than the winter value, strict comparison 
between the summer and winter k-values cannot be made for two reasons. 
Firstly, the point of tracer release was different for summer and 
winter. Secondly, the lateral mixing process depends strongly on bend 
induced spiral motions and the intensity of these motions together 
with the size of the dominant turbulent eddies vary laterally. Conse- 
quently tracer experiences different rates of mixing depending on the 
source position. Another factor which should be mentioned is that, in 
the laboratory, comparison of Ew or k values from ice-cover and 
open water tests were made for corresponding flows having approximately 


the same depth and discharge but different slopes. The summer test 


was performed at a higher depth and discharge compared to the winter 
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test although the slope during both tests is not significantly different. 


In spite of the above limitations, some rough comparisons can 
be made between the winter and summer results. The mixing potential of 
the river as measured by Ee. for the initial reach length of about 
7000 ft is more than four times greater in summer than in winter. This 
can be seen by comparing the concentration distributions (Figures 7.5a 
and b) at Sections 3 (winter) and E (summer) which are located closely 
together and indicate a similar degree of mixing, even though in summer 
the tracer was injected near the left bank, while injection was near 
the channel centerline in winter. The low value of Ew for winter is 
reasonable and is probably caused by the reduction in both diffusion 
capacity and strength of spiral motions as discussed previously in 


Chapter 6. 


It is important to remember that E ow is considered to be 
representative of the reach average of the generally variable os 
(depth-averaged transverse mixing coefficient). It was emphasized in 
Chapter 3 that if >. (or Es does vary in a particular flow situa- 
tion, different values of Ey. or k will be obtained depending on 
the point (2) about which the moments are taken. It is to be noted 
that both laboratory and field data were analysed with Z 5 taken as 
the transverse co-ordinate of source position. It is of particular 
interest, therefore, to test the sensitivity of k to moments taken 
about different Ze Table 7-4 summarises the variation of the 


normalised exchange ocefficient, k , for x 0.0 and - 6.0. 
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Also shown in Table 7-4 are the k-values obtained by taking moments 
about the centroid of the tracer flux distribution as suggested by 


Holley and Abrahams (1972) for centerline release of tracer. 


It is important to mention that the winter test achieved 90% 
mixing over the test reach with the tracer cloud occuping the entire 
river cross-section. Table 7.4 indicates significant differences in 
the computed values of k for the various reaches but the distance 
weighted average k for the entire reach resulting from taking moments 
about z=0 and z=- 6.0 are practically the same. The reason 
for this is probably due to the fact that the cumulative effects of 
variations in oe (or o) which were not accounted for in Equation 7.7 
average out when the mixing is considered over a series of alternating 


bends and when a very high degree of mixing has been achieved. 


Table 7-4 also indicates that approximately the same average 
k-value is obtained when moments are taken about the centroid of the 
tracer flux distribution. This seems fortuitous because the locus of 
the centroid of the tracer flux distribution at the various sections 
is neither concentric with or parallel to the channel axis which is 


taken either as an arc of a circle or a straight line. 


Figure 7.8 indicates that the longitudinal variation of the 
normalised transverse mixing coefficient for the winter test follows 
a trend which is quite similar to that observed in the laboratory 


tests with the maximum k occurring around the middle of a bend and 
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TABLE 7.4 DEPENDENCE OF E og! Ue ON POINT ABOUT WHICH 
MOMENTS OF TRACER FLUX ARE TAKEN OR VALUE OF 
CUMULATIVE DISCHARGE ABOUT WHICH MOMENTS OF 
TRACER CONCENTRATION ARE TAKEN. 


Cross- 
Section 


Number. 
Moments Moments 


About About 


Distance 
Weighted 
Averaged, 


NOTES. 


+ E og! Ue obtained from moments of tracer 


distribution with respect to gz. 
++ E og! tx obtained from moments of tracer 
concentration distribution with respect 
* 2 is transverse coordinate of centroid 


tracer flux distribution. 
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FIGURE 7.8 LONGITUDINAL VARIATION OF THE NORMALISED 
TRANSVERSE MIXING COEFFICIENT, Eg! tak. 
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the minimum at the beginning of a bend. 


The open water tests indicated some reduction in the value 
of normalised mixing coefficient, k , between sections C and D. This 
could be expected again because the transverse velocities would tend 
to move the tracer back toward the left bank and thus inhibit the 
spreading of tracer which was injected from the left bank. It is noted 
that the open water results are not as reliable as those for the winter 


test. 


Yotsukura and Cobb's (1972) method for describing the tracer 
concentration distribution as a function of the normalised cumulative 
discharge, that is, c= e ¢q,/0), appears to be more suitable for non- 
uniform meandering channels such as the Lesser Slave River. It is to 


be noted that q, is the cumulative discharge and Q is the total 


discharge. The cumulative discharge, q, , is defined as: 


where z= 0 is taken at one bank. 


Representing the concentration distribution as c= © (q,/Q) 
instead of c= c (z/w) describes the mixing process in the context 
of fluid mixing rather than spatial transport. By using q, as an 
alternative transverse coordinate in place of z , Yotsukura and Cobb 
showed that the steady state convective-diffusion equation averaged 


over depth could be transformed to: 
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It is noted that o. = dcx, q,) and as mentioned previously, it 
defines the variation of the depth-averaged transverse mixing coef- 


ficient for a given flow situation. 


The moments technique can now be applied to Equation 7.9, 
again assuming that o, = 1.0, by taking moments with respect to 4, ’ 
where q, cai! a qd, and qd. is the cumulative discharge about which 
the moments are taken. The resulting equation after some simplifica- 


tions is: 


q q = 
2 , 2 mn 
eq. d C¢ —_~ dq 
“do SI 94, ) 
q q 
- eee = 2E ow eA a Pe (7.10) 
oa a q2 
| cdq, | cdq) 
qy qy 


in which qy and q, are the cumulative discharges measured from q,70- 
Equation 7.10 can be rewritten as: 


aM, 


dX 


E ow Bo (7.11) 
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where the definitions of a and BLO) can be obtained by compa- 


ring Equations 7.11 and 7.10. 


Equation 7.11 was applied to the winter data and the moments 
were taken about the point of tracer release, that is Sptle 0.44 Q. 
The numerical values of Nye eS é Che) and ade are given in Table 
7-5. The overall average transverse mixing coefficient was evaluated 
as 0.1085 and the corresponding E og! tak as 0.34. It is seen that the 
overall average values of 2 pao obtained from Equation 7.7 by correct- 


ing for the moments and Equation 7.11 are not significantly different. 


TABLE 7-5 NUMERICAL VALUES OF My)» PCS, AND Eee IN EQUATION 
Perl 


oes 


Distance 


Section aaere 


-2 
feo aa 


21882 
39208 
60689 
68886 
83201 


Average 


Again, because a5 varies both longitudinally and laterally 
different values of Ew are obtained depending on the value of A. 


about which the moments are taken. Table 7.4 gives a comparison be- 
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tween the values of E oy! Ua® evaluated by taking moments about ae 


0.44 Q and aoe 0.55 Q; 


It is evident from the above discussions that if the mixing 
is considered over a series of bends and provided a high degree of 
mixing has been achieved then the method of moments (Equation 7.7) 
combined with the adjustment (Equation 7.8) proposed by Fischer give 
satisfactory estimates of the overall k irrespective of the point 


about which the moments are taken. 


Flow in a channel of non-uniform cross-section is usually 
characterised by an average velocéity and depth. For this situation 
the transverse mixing coefficient can be obtained by applying the 


formula: 


(7412) 


+ ad . . i - ; 
where Oo, is the variance of the concentration distribution. This 


equation, as outlined by Sayre and Chang (1968), is based on the assump- 


tions of a gradient-type mixing and a uniform downstream velocity at 
each point on the channel cross-section. Equation 7.12 is valid only 
until a significant amount of tracer encounters the opposite bank for 


bank release or both banks for centerline release of tracer. 


When Equation 7.12 was used to determine the mixing coef- 
ficient for the summer test, it resulted in an average value of Ew = 


0.86 ee fase and k = 0.65. This k-value is almost 100% larger than 
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the value predicted by Equation 7.7 which takes account of depth and 
velocity variations in the lateral direction. That the assumptions of 
constant depth and uniform velocity give unsatisfactory predictions of 
k is also confirmed by Holley et al (1972). They demonstrated numeri- 
cally that if a trapezoidal channel is represented by a rectangular 
channel having the same area and surface width, the mixing coefficient 


calculated from the concentration distributions may be larger. 


Reports of field measurements of the normalised transverse 
mixing coefficient are listed in Table 7.6. It is evident from this 
table that the k-values obtained are within the same range as those 


observed for idealised meandering laboratory channels. 


The results of Holley and Abrahams (1973) particularly de- 
monstrate how increased mixing caused by obstructions can increase the 
magnitude of k. Their test reaches had relatively mild bends but 
they were bordered by groins which, according to their laboratory 


study, can, increase: & “by a@ Lactor tof 2 to 3. 


It was shown by dimensional arguments in Section 6.2.4 that: 
E f 

L 
satchel T z. = a | (6027) 


An indication of the importance of r/R can be seen in the two 
results obtained for the Missouri River. It is assumed that d, is 


approximately equal to R for open water conditions. 
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TABLE 7.6 SOME FIELD MEASUREMENTS OF THE NORMALISEDTRANSVERSE . 


SOURCE 


GLOVER (1964) 


YOTSUKURA et al 
(1970) 


FISCHER (1967) 


YOTSUKURA and 
COBB (1972). 


HOLLEY and 
ABRAHAM (197 3) 


SAYRE and YEH 
(1973) 


HIGHWAY and RIV- 
ER ENGINEERING 
DIVISION ALBERTA 
RESEARCH (1974) 


PRESENT INVESTI- 
GATION 


* 
This Result Is A Preliminary Estimate. 


MIXING COEFFICIENT. 


Columbia River 
Near Richland 


Missouri River 
Near Blair, 
Nebraska. 


Atrisco Feeder 
Canal. 


South River. 


Bernando Convey- 
ance Canal 


Ijssel River. 


Waal River. 


Missouri River 
Near Brownville 
Nebraska, 


Athabasca River 
Downstream of 
Fort McMurray 
( ICE-COVERED) 


Lesser Slave 
River. 


Lesser Slave 


River (ICE-COVERED 


Test Reach with 
a Bend of Radius 
12000 ft. 


Reasonably 
Straight Reach. 


Test Reach had a 
few Mild Bends. 


Straight Reach 
with Shifting 
Bed Configuration 


River Bordered by 
Groins. 


River Bordered by 
Groins. 


Test Reach had 2 
Bends with Radii 
6400 and 3400 ft. 


ild Bends with 
Some Islands. 
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The curvature of the Brownville reach was about 2 to 3 more pronounced 
than that of the Blair reach with the r/R ratio 3 to 5 times greater 
in the Brownville test. According to Sayre and Yeh (1973) the river 
discharge at the time of the Blair experiment was approximately 40% 
less than that during the Brownville test. However, they attributed 
the greater rate of transverse mixing in the Brownville test more to 
difference in sinuosity of the two reaches than to difference in river 
discharge. This seems reasonable since the friction factor, f , for 
both reaches was adunox tiara the same and the aspect ratio was not 


significantly different. 


The effect of width to depth ratio is not clear from Table 
7-6. However, if data obtained for rivers with mild bends and no 
obstructions are compared then the data seem to indicate some depen- 


dence of k on the aspect ratio, with k increasing with W/R. 
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CHAPTER 8 


SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


Transverse diffusion and mixing in straight and meandering 
ice-covered and open channel flows were investigated by means of a non- 
buoyant tracer released into the flow at a steady and continuous rate. 
Tracer concentration and velocity measurements were made in the labora- 
tory for similar conditions of discharge and depth in both open and 
ice-covered channels. Transverse mixing experiments for both open and 
ice-covered conditions were also conducted iy the Lesser Slave River. 
The laboratory and field results were analysed and used to evaluate 
the effect of an ice-cover on the transverse diffusion and mixing co- 


efficients. 


8.1 Conclusions 


The primary aim of this thesis was to demonstrate the signi- 
ficant reduction in the mixing capacity of a channel (as measured by 
E. or Ew caused by an ice-cover. The laboratory and field tests 
indicate that this reduction in absolute terms can vary from 50 to 500%. 
However, when the mixing coefficients are normalised ae the product 


u,R reductions of up to 300% could be observed. Other findings of 


this investigation are as follows. 


8.2 Results Related to Straight Flume 


1. The depth-averaged transverse concentration profiles at various 
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downstream stations for a given flow condition can be represented 
by a single curve if Cees is plotted against z/b. This curve 


was Gaussian for both open and ice-covered tests. 


From the above suggestion of similarity, a formula was derived for 
evaluating the transverse diffusion coefficient. This formula, 


2 
E. = USP /4xln2 , requires only one measured depth-averaged trans- 


verse concentration profile and was found to give satisfactory 


results. It is only applicable to situations with the source 


located at the channel centerline. 


For side release of tracer, the transverse diffusion coefficient 
was found to increase with distance from the source and approached 
the constant value for centerline injection as the tracer cloud 


occupied an increasing area of the uniform flow region. 


For similar conditions of discharge and depth the presence of an 
ice-cover was observed to reduce the diffusion capacity, as mea- 


sured by E. » by as much as one-half the value for open water 


conditions. 


The values of the normalised transverse diffusion coefficient, k 
(= E,/u,R), for ice-cover conditions were generally found to be 

slightly less than those for the corresponding open channel tests. 
These latter k-values are within the range of previously reported 


k-values in open channel flows. 
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For fully rough turbulent flow k was found to increase with in- 
crease in f£ when W/R was held constant. No obvious dependence 


of k on W/R was indicated by the available data. 


Results Related to Meandering Flume 


The growth of the second mements of the tracer flux distribution 


-with longitudinal distance was non-linear. The growth rate of 


M, was observed to decrease and in some situations became nega- 


2 


tive, particularly where the spiral motion decayed and reversed. 


The local depth-averaged mixing coefficient, ree » was found to 
vary both laterally and longitudinally and seemed to be associated 


largely with variations in the intensity of the spiral motions. 


The magnitude of the averaged normalised mixing coefficient, k , 
over the test reach for ice-covered flows was observed, on the 
average, to be about 2 times less than the corresponding value in 
open channel flows for similar flow conditions. The range of 
k-values obtained are comparable to those reported in the liter- 


ature for field conditions. 


The magnitude of the transverse mixing coefficient, E ow » was 
observed to be of the order 4 times smaller for ice-cover condi- A 
tions than for the corresponding open channel. This substantial 


reduction is presumably due to reduction in the scale of the domi- 
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nant turbulent eddies and strength of the spiral motions in the 


presence of an ice-cover. 


An apparent dependence of the transverse mixing coefficient on 
source position was observed. It was suggested that it is pro- 
bably a manifestation of the effects of transverse velocities on 


the mixing. The mixing coefficient varied longitudinally for all 


source positions and the averaged lateral mixing coefficient for 


different source positions was found, in general, to reach a 
maximum within the middle third of a bend and a minimum at the 
beginning of a bend. Again this variation would seem to be related 


to the growth, decay and reversal of the spiral motions. 


Negative values of the transverse mixing coefficient were sometimes 
observed, particularly at the beginning of a bend where there is 
reversal of the spiral motion. A negative ae value implies 
convective transport by spiral motions of low depth-averaged con- 


centration fluid to regions of high depth-averaged concentration. 


The principal parameters controlling the normalised mixing coef- 
ficient, k , were identified as the aspect ratio W/d, (or W/R), 
r/d, and: they frictton factor, ‘£. For the present laboratory 
investigation, k was found to decrease with increasing r/d, 


ratio, increasing roughness and decreasing aspect ratio. 


Of the two methods proposed for evaluating the transverse mixing 


coefficient, the integral method was found to give extremely high 
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k-values for situations where the tracer was released at the chan- 
nel centerline. The cause of this could not be clearly identified. 
Both integral and moments methods yielded similar patterns for the 
longitudinal variation of E ow’ The overall average Bae for a 
given source location as computed by the integral method was on 
the average only about 25% higher than the corresponding value 
estimated by the moments method. However, large differences in 
_the value of Ee for individual sub-reaches were sometimes ob- 


served between the two methods. 


8.4 Results Related to Field Investigations 


Since the present data relate only to field measurements on 
one stream, general conclusions would be premature but the results do 
indicate that an ice-cover can reduce the mixing capacity of a river 
substantially. This is confirmed by the laboratory data. The assump- 
Eon of quick vertical mixing normally made in the context of field 
studies was found not to be valid in the case of relatively deep, 


meandering rivers. This was also found to be true in the laboratory 


tests. 


8.5 Recommendations 


No systematic and detailed investigations have so far been 


carried out to delineate the effects of various parameters controlling 
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the transverse spreading of tracer, even for the idealized conditions 


in the laboratory. The following recommendations should be considered 


in planning future research on transverse mixing and diffusion: 


LBA 


Detailed measurements of the normalised mixing and diffusion 
coefficients along lines similar to that adopted by Miller 

and Richardson (1974) are needed. These measurements should 
be used to produce a Moody type diagram which should complete- 
ly describe the dependence of E_/u,R or E yg! Ua® on other 


parameters such as W/R, f and r/R. 


Present knowledge of the hydraulics of flow in ice-covered 
channels is very limited. Since the mixing phenomena are 
closely related to the velocity field, detailed velocity mea- 
surements, particularly in ice-covered meandering channels 
are necessary for a better understanding of mixing in channel 


flows. 
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APPENDIX A 
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FIGURE A.1lb VERTICAL DISTRIBUTIONS OF HORIZONTAL ANGLE OF 
DEVIATION OF VELOCITY VECTOR, Runs 401/2/3, 
SECTIONS: seawae Si/, TO. By. 1/2. 
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FIGURE A.llc VERTICAL DISTRIBUTIONS OF HORIZONTAL ANGLE OF 
DEVIATION OF VELOCITY VECTOR, Runs 401/2/3, 
SECTIONS 8, = Wa To 5/6. 
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FIGURE A.12a VERTICAL DISTRIBUTIONS OF HORIZONTAL ANGLE OF 
DEVIATION OF VELOCITY VECTOR, Runs 407/8, 
SECTIONS fee "ig TO S, 2 Aaw/s. 
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VERTICAL DISTRIBUTIONS OF HORIZONTAL ANGLE OF 
DEVIATLON OF VELOCITY VECTOR, Runs 407/3, 
SECTIONS 9, = Su/g To 9, = Fz 


ee 


oe - 
f k a 
i ; 
z — ant = So. = 


! er | 
a  @ + Sea G ° 


di : aobigeb ni Y 7 


| ei | 0 


a : re ote r 
Te Spas) ¥ j ny 
oA cay Y , 


: : 7 ee a . 2 - a ey 
ra 4d . ’ : 9 = aa nd ait ries Dee 


io - ee & 1 < iy 
oe i? wee 
: A P ‘ak 
rH =) Sa &  “~ Ogg re 
; oe aly ae 
Ne 7 rs " 


S\ 4 ; Le ne 
Ta ee ee 
‘a . Pad a @ ; ” yi ae a 
7 == ai Me : » ; : - ae 4 : a : a4 = ' ! i * ’ of "5 
; 7 3 a 1 i. 7 a he 7 : > =T a Baie 
wae ~~ -) rr en eect) ee ea 


i eh 
a — Vea) ei) & i) r Q 
Aiayiget a IE 
Wr atin Was 
or yes i} ae F 
hed en oe q Pe fy 
a ene 
dg, 
if 
iF \ 
i 7 
me | on 
& j or i 
| ees 
abs hs 
i, oh 
a fir ae Mit 
a i ae 
r re be 
ae t j 
oeer Cie 
\ - zi 
y ig ge iW eo 
| | } 
j 


& Wy ; is +h 
a 


pees 
ne 


GOR crit alata al 


= 


ae > ru a 
7 : a é «| BN Wu 
7 Par 


“a2 
<a Pi 


D 
Lee) 


a pee 
- a 
— <a. - 
S 
—— : 


